System and method for acoustic echo cancelation using deep multitask recurrent neural networks

ABSTRACT

A system for performing echo cancellation includes: a processor configured to: receive a far-end signal; record a microphone signal including: a near-end signal; and an echo signal corresponding to the far-end signal; extract far-end features from the far-end signal; extract microphone features from the microphone signal; compute estimated near-end features by supplying the microphone features and the far-end features to an acoustic echo cancellation module including a recurrent neural network including: an encoder including a plurality of gated recurrent units; and a decoder including a plurality of gated recurrent units; compute an estimated near-end signal from the estimated near-end features; and transmit the estimated near-end signal to the far-end device. The recurrent neural network may include a contextual attention module; and the recurrent neural network may take, as input, a plurality of error features computed based on the far-end features, the microphone features, and acoustic path parameters.

CROSS-REFERENCE TO RELATED APPLICATION(S)

This application is a continuation-in-part of U.S. patent application Ser. No. 16/573,573, filed in the United States Patent and Trademark Office on Sep. 17, 2019, which claims priority to and the benefit of U.S. Provisional Patent Application No. 62/825,681, filed in the United States Patent and Trademark Office on Mar. 28, 2019 and U.S. Provisional Patent Application No. 62/838,146, filed in the United States Patent and Trademark Office on Apr. 24, 2019, and this application claims priority to and the benefit of U.S. Provisional Patent Application No. 62/914,875, filed in the United States Patent and Trademark Office on Oct. 14, 2019, the entire disclosures of which are incorporated by reference herein.

FIELD

Aspects of embodiments of the present disclosure relate to the field of audio signal processing, including specifically a system and method for providing acoustic echo cancellation using deep multitask recurrent neural networks.

BACKGROUND

During a real-time or near real-time audio or videoconferencing interactions, a far-end participant may notice acoustic echoes in the signal received from a near-end participant. Acoustic echoes may be generated when a far-end signal playing out of a near-end loudspeaker is coupled back to a microphone at the near-end. When the signal from the near-end microphone is transmitted to the far-end, the far-end user hears a mixture of a near-end signal (e.g., speech from the near-end participant) and a delayed and modified version of their own far-end voice as the acoustic echo. An acoustic echo canceller (AEC) or acoustic echo suppressor (AES) attempts to reduce or remove this echo while avoiding distorting the speech of the near-end user (e.g., without distorting the near-end signal).

SUMMARY

Aspects of embodiments of the present disclosure relate to cancelling acoustic echoes in an audio communication system using deep multitask recurrent neural networks.

According to one embodiment of the present disclosure, a system for performing echo cancellation includes: a processor; and memory storing instructions that, when executed by the processor, cause the processor to: receive a far-end signal from a far-end device at a near-end device; record a microphone signal at the near-end device, the microphone signal including: a near-end signal; and an echo signal corresponding to the far-end signal; extract a plurality of far-end features from the far-end signal; extract a plurality of microphone features from the microphone signal; compute a plurality of estimated near-end features from the microphone signal by supplying the microphone features and the far-end features to an acoustic echo cancellation module including a contextual attention neural network including: an encoder including a plurality of gated recurrent units configured to encode the microphone features and the far-end features as encoded features; a contextual attention module configured to compute contextual attention features from the encoded features; and a decoder including a plurality of gated recurrent units configured to decode the estimated near-end features from the contextual attention features; compute an estimated near-end signal from the estimated near-end features; and transmit the estimated near-end signal to the far-end device.

The encoder may include one or more of layers of gated recurrent units, and the encoded features may correspond to outputs of a last layer of the one or more layers of the gated recurrent units of the encoder, the encoded features may be supplied to a multi-head self-attention layer of the contextual attention module; and the decoder may include one or more layers of gated recurrent units, and the estimated near-end features correspond to an output of a last layer of the one or more layers of gated recurrent units of the decoder.

The contextual attention module may further include a multi-head attention layer configured to receive an output of the multi-head self-attention layer.

The contextual attention module may further include a first layer normalization layer between the multi-head self-attention layer and the multi-head attention layer.

The contextual attention module may further include a second layer normalization layer configured to compute the contextual attention features in accordance with outputs from other layers of the contextual attention module.

The memory may further store instructions that, when executed by the processor, cause the processor to: compute a plurality of acoustic path parameters in accordance with prior acoustic path parameters, prior far-end features, and prior microphone features; and compute a plurality of error features based on the far-end features, the microphone features, and the acoustic path parameters, and the acoustic echo cancellation module may further be configured to compute the plurality of estimated near-end features from the microphone signal based on the error features.

The far-end features supplied to the contextual attention neural network may include a current frame of far-end features and a causal window of a plurality of previous frames of far-end features, and the microphone features supplied to the contextual attention neural network may include a current frame of microphone features and the causal window of a plurality of previous frames of microphone features.

The estimated near-end features may include a current frame of estimated near-end features and the causal window of a plurality of previous frames of estimated near-end features, and the contextual attention neural network may be trained by iteratively training a plurality of parameters configuring contextual attention neural network to minimize differences between the estimated near-end features and a plurality of ground truth near-end features of the plurality of training data.

The training data may be generated by: loading a corpus of training data including recorded utterances from a plurality of different human speakers; selecting a plurality of pairs of human speakers; for each pair of human speakers: concatenating multiple utterances of a first speaker of the pair of human speakers to generate a training far-end signal; transforming the concatenated utterances to simulate an acoustic path to generate a training echo signal; padding an utterance of a second speaker of the pair of human speakers to generate a training near-end signal having a length equal to the training far-end signal; and mixing the training echo signal with the training near-end signal to generate a training microphone signal.

The transforming the concatenated utterances may include convolving the far-end signal with a room impulse response of a simulated room.

The transforming the concatenated utterances may further include: applying hard clipping to the far-end signal to generate a clipped far-end signal; and applying sigmoidal distortion to the clipped far-end signal.

The far-end features, the microphone features, and the estimated near-end features may include log short time Fourier transform features in logarithmic spectral space.

According to one embodiment of the present disclosure, a system for performing echo cancellation includes: a processor; and memory storing instructions that, when executed by the processor, cause the processor to: receive a far-end signal from a far-end device; record a microphone signal including: a near-end signal; and an echo signal corresponding to the far-end signal; extract a plurality of far-end features from the far-end signal; extract a plurality of microphone features from the microphone signal; compute a plurality of acoustic path parameters in accordance with prior acoustic path parameters, prior far-end features, and prior microphone features; compute a plurality of error features based on the far-end features, the microphone features, and the acoustic path parameters; compute a plurality of estimated near-end features from the microphone signal by supplying the microphone features, the far-end features, and the error features to an acoustic echo cancellation module including a recurrent neural network including: an encoder including a plurality of gated recurrent units configured to encode the microphone features and the far-end features as encoded features; and a decoder including a plurality of gated recurrent units configured to decode the estimated near-end features based on the encoded features; compute an estimated near-end signal from the estimated near-end features; and transmit the estimated near-end signal to the far-end device.

The acoustic echo cancellation module may further include a contextual attention module configured to compute contextual attention features from the encoded features, and the decoder may be configured to compute the estimated near-end features based on the contextual attention features.

The far-end features supplied to the acoustic echo cancellation module may include a current frame of far-end features and a causal window of a plurality of previous frames of far-end features, the microphone features supplied to the acoustic echo cancellation module may include a current frame of microphone features and the causal window of a plurality of previous frames of microphone features, and the error features supplied to the acoustic echo cancellation module may include a current frame of error features and the causal window of a plurality of previous frames of error features.

The estimated near-end features may include a current frame of estimated near-end features and a causal window of a plurality of previous frames of estimated near-end features, and the recurrent neural network may be trained by iteratively training a plurality of parameters configuring the recurrent neural network to minimize differences between the estimated near-end features and a plurality of ground truth near-end features of the plurality of training data.

The training data may be generated by: loading a corpus of training data including recorded utterances from a plurality of different human speakers; selecting a plurality of pairs of human speakers; for each pair of human speakers: concatenating multiple utterances of a first speaker of the pair of human speakers to generate a training far-end signal; transforming the concatenated utterances to simulate an acoustic path to generate a training echo signal; padding an utterance of a second speaker of the pair of human speakers to generate a training near-end signal having a length equal to the training far-end signal; and mixing the training echo signal with the training near-end signal to generate a training microphone signal.

The transforming the concatenated utterances may include convolving the far-end signal with a room impulse response of a simulated room.

The transforming the concatenated utterances may include: applying hard clipping to the far-end signal to generate a clipped far-end signal; and applying sigmoidal distortion to the clipped far-end signal.

The far-end features, the microphone features, and the estimated near-end features may include log short time Fourier transform features in logarithmic spectral space.

BRIEF DESCRIPTION OF THE DRAWINGS

The accompanying drawings, together with the specification, illustrate exemplary embodiments of the present disclosure, and, together with the description, serve to explain the principles of the present disclosure.

FIG. 1A is a schematic diagram illustrating a communication system including an acoustic echo cancellation (AEC) module configured to reduce or remove acoustic echoes.

FIG. 1B is a schematic diagram illustrating a communication system including an acoustic echo cancellation (AEC) module and a residual echo suppression (RES) module configured to reduce or remove acoustic echoes.

FIG. 2 is a block diagram illustrating an acoustic echo cancellation system, according to one embodiment of the present disclosure.

FIG. 3 is a block diagram illustrating an input layer of an unrolled deep gated recurrent unit (GRU) network, according to one embodiment of the present disclosure.

FIG. 4 is a block diagram illustrating an architecture of an acoustic echo cancellation (AEC) neural network, according to one embodiment of the present disclosure.

FIG. 5A is a flowchart depicting a method for training a recurrent neural network of an acoustic echo cancellation module, according to one embodiment of the present disclosure.

FIG. 5B is a flowchart of a method for generating training data, according to one embodiment of the present disclosure.

FIG. 5C is a flowchart of a method for generating a simulated echo signal y(t), according to one embodiment of the present disclosure.

FIG. 6A is a block diagram illustrating an architecture of an acoustic echo cancellation (AEC) neural network operating in inference mode, according to one embodiment of the present disclosure.

FIG. 6B is a flowchart depicting a method for computing an estimated near-end signal from a received far-end signal and a microphone signal using a neural network in inference mode, according to one embodiment of the present disclosure.

FIG. 7 is a block diagram illustrating a hybrid system including the deep multi-task acoustic echo cancellation system, according to one embodiment of the present disclosure.

FIG. 8A is a schematic block diagram of an acoustic echo cancellation system in which a deep neural network (denoted as “EchoDNN”) is used to cancel echoes from a microphone signal, according to one embodiment of the present disclosure.

FIG. 8B is an architectural block diagram of a deep neural network configured to cancel echoes from a microphone signal, according to one embodiment of the present disclosure.

FIG. 9 is a block diagram illustrating an acoustic echo cancellation system using a generative network, according to one embodiment of the present disclosure.

FIG. 10 is a block diagram illustrating an acoustic echo cancellation system including a contextual attention neural network, according to one embodiment of the present disclosure.

FIG. 11A is a block diagram depicting a contextual attention aware neural network, according to one embodiment of the present disclosure.

FIG. 11B is a more detailed block diagram depicting layers of a contextual attention neural network, according to one embodiment of the present disclosure.

FIG. 11C is a block diagram depicting a contextual attention module, according to one embodiment of the present disclosure.

FIG. 12 is a flowchart depicting a method for computing an estimated near-end signal from a received far-end signal and a microphone signal using a contextual attention neural network in inference mode, according to one embodiment of the present disclosure.

FIG. 13 is a block diagram illustrating an acoustic echo cancellation system including a hybrid of a Frequency Domain Normalized Mean Square (FDNLMS) technique in conjunction with a contextual attention neural network, according to one embodiment of the present disclosure.

FIG. 14A is a block diagram depicting a contextual attention neural network, where the contextual attention neural network is configured to accept far-end features, error features, and microphone features as input, according to one embodiment of the present disclosure.

FIG. 14B is a more detailed block diagram depicting layers of a contextual attention neural network configured to accept far-end features, error features, and microphone features as input, according to one embodiment of the present disclosure.

FIG. 15 is a flowchart depicting a method for computing an estimated near-end signal from a received far-end signal and a microphone signal using a hybrid of a classical acoustic echo cancellation technique and a contextual attention neural network in inference mode, according to one embodiment of the present disclosure.

FIG. 16 is a block diagram of a near-end device, according to one embodiment of the present disclosure.

DETAILED DESCRIPTION

In the following detailed description, only certain exemplary embodiments of the present disclosure are shown and described, by way of illustration. As those skilled in the art would recognize, the disclosure may be embodied in many different forms and should not be construed as being limited to the embodiments set forth herein. In the figures and the below discussion, like reference numerals refer to like components.

FIG. 1A is a schematic diagram illustrating a communication system including an acoustic echo cancellation (AEC) module configured to reduce or remove acoustic echoes. As shown in FIG. 1A, a far-end signal x(t) is received and played back on a near-end loudspeaker 12 of a near-end device 10. Generally, the far-end signal x(t) contains speech from a far-end user or participant in the interaction, where the far-end signal x(t) is captured by a far-end communication device which transmits the far-end signal x(t) to be received by the near-end communication device (e.g., over a communications network such as a cellular communications network). As such, for the sake of convenience, the far-end signal x(t) may sometimes be referred to herein as “far-end speech.” The sound produced by the near-end loudspeaker 12 may be detected by a near-end microphone 14 as acoustic echoes y(t), and the sound may travel to the near-end microphone 14 on multiple different echo paths 16, such as after the sounds bounce off walls and other surfaces 18 in a physical environment in which the near-end device 10 is located. The near-end microphone 14 may also detect desired near-end signal s(t), such as speech from a near-end user or participant in the interaction. The signal d(t) produced by the near-end microphone 14 therefore may contain both the desired near-end signal s(t) and the undesired acoustic echoes y(t). The microphone signal d(t) is provided to the acoustic echo cancellation (AEC) module 100, which attempts to remove the echo signal y(t) and to produce an output signal q(t) that approximates the desired near-end signal s(t) (e.g., corresponds to the microphone signal d(t) with the acoustic echoes y(t) removed). The processed the output signal q(t) can then be sent to the far-end as an estimate of the desired near-end signal s(t) with the echo signal y(t) removed.

For the sake of clarity, as used herein, given v(t) as an arbitrary time-domain signal at time t: the short-time Fourier transform (STFT) complex-valued spectrum of v(t) at frame k and frequency bin f is denoted by V_(k,f); its phase is denoted by ∠V_(k,f); and its logarithmic magnitude is denoted by {tilde over (V)}_(k,f). {tilde over (V)}_(k) represents the vector of logarithmic magnitudes at all frequency bins f and frame k.

FIG. 1B is a schematic diagram illustrating a communication system including an acoustic echo cancellation (AEC) module and a residual echo suppression (RES) module configured to reduce or remove acoustic echoes. The microphone signal d(t) includes near-end speech signal s(t) and acoustic echo y(t): d(t)=s(t)+y(t) In some embodiments, the microphone signal d(t) also includes other components such as additive noise n(t) (e.g., d(t)=s(t)+y(t)+n(t)). The acoustic echo signal y(t) is a modified version of far-end speech signal x(t) and includes room impulse response (RIR) and loudspeaker distortion, both of which may cause nonlinearities in the relationship between x(t) and y(t).

Broadly, the acoustic echo cancellation (AEC) problem is to retrieve the clean near-end signal s(t) after removing acoustic echoes due to detection of the far-end signal x(t) by the near-end microphone 14. Comparative systems, as shown in FIG. 1B, estimate a model of the echo path with a linear adaptive filter 110 and then subtract the estimated echo v(t) from the microphone signal d(t). In addition, to enable echo estimation, some comparative AEC methods use a double-talk detector (DTD) to halt the filter adaptation during double-talk periods, when both the near-end and far-end signals are simultaneously present (in contrast to single-talk periods where only the near-end or only the far-end signal is present or non-negligible). Usually, after applying the AEC, there still exists some amount of residual echo at the output of the AEC method (d(t)−v(t)). One of the reasons that comparative AEC methods cannot perfectly remove all the echo noises is that the echo path 16 is not a linear function of the far-end speech signal x(t), even when the echo path 16 is perfectly linear. For example, the audio amplifier (e.g., a power amplifier) and near-end loudspeaker 12 can cause nonlinearities in the relationship between the far-end speech signal x(t) and the echo signal y(t). In addition, the length of the adaptive linear filter may not be sufficient to remove longer echoes (e.g., in large rooms or hallways). Accordingly, a residual echo suppressor (RES) 150 can be further applied to improve the near-end signal, where the RES may be realized by Wiener filter or spectral subtraction in the frequency domain. The final output of AEC system is estimated near-end signal q(t).

Aspects of embodiments of the present disclosure relate to the recurrent neural network (RNN) architectures for acoustic echo cancellation (AEC). Some embodiments relate to the use of deep gated recurrent unit (GRU) networks (see, e.g., K. Cho, B. van Merrienboer, C. Gulcehre, D. Bahdanau, F. Bougares, H. Schwen, and Y. Bengio, “Learning Phrase Representations using RNN Encoder-Decoder for Statistical Machine Translation,” in Proc. Empirical Methods in Natural Language Processing, 2014, pp. 1724-1734. and J. Chung, C. Gulcehre, K. Cho, and Y. Bengio, “Empirical evaluation of gated recurrent neural networks on sequence modeling,” in Proc. NIPS Deep Learning Workshop, 2014) in an encoder/decoder architecture to map the spectral features of the microphone signals d(t) and far-end signals x(t) to a hyperspace (e.g., a feature space such as logarithmic spectral space), and then decode the target spectral features of the near-end signal s(t) from the encoded hyperspace. In some embodiments, the RNN acoustic echo cancellation module is trained using multitask learning to learn an auxiliary task of estimating the echo signal y(t) in order to improve the main task of estimating the clean near-end speech signal s(t) as estimated near-end signal q(t). As discussed in more detail below, experimental results show that embodiments of the present disclosure cancel acoustic echo in both single-talk and double-talk periods with nonlinear distortions without requiring a separate double-talk detector.

FIG. 2 is a block diagram illustrating an acoustic echo cancellation system, according to one embodiment of the present disclosure. A microphone signal d(t)=s(t)+y(t) is supplied to an acoustic echo cancellation system 200 to cancel or remove the echo signal y(t). For the sake of convenience, aspects of embodiments of the present disclosure will be described herein in the context of audio signals that are sampled at 16 kHz. However, embodiments of the present disclosure are not limited thereto, and may be applied with audio signals sampled at higher rates than 16 kHz or at lower rates than 16 kHz.

In the embodiment shown in FIG. 2, the microphone signal d(t) is supplied to a microphone signal feature extraction module 210, which includes a short time Fourier transform (STFT) module 212, an absolute value module 214 to remove the conjugate symmetric half of the output of the STFT module 212, and a logarithmic operation module 216 to compute the final logarithmic magnitude spectral feature vector or logarithmic spectral features or microphone signal features {tilde over (D)}_(k,f) in “feature space” or “hyperspace” or “logarithmic spectral space.” The STFT module 212 may also compute a phase ∠D_(k,f) of the microphone signal d(t).

For the sake of convenience, aspects of embodiments of the present disclosure will be described herein where the spectral feature vectors are computed using a 512-point short time Fourier transform (STFT) with a frame shift of 256-point (given the 16 kHz sampling rate, each frame corresponds to 32 milliseconds with a 16 millisecond shift between frames, resulting is 16 milliseconds of overlap between frames). In some embodiments, the absolute value module 214 reduces the 512-point STFT magnitude vector to 257-point by removing the conjugate symmetric half. In some embodiments, the features (e.g., the microphone signal features {tilde over (D)}_(k,f)) are standardized to have zero mean and unit variance using the scalars calculated from the training data, as discussed in more detail below. As would be understood by one of skill in the art, the spectral feature vectors may be computed with more than 512 points or fewer than 512 points and with longer or shorter frame shifts (e.g., more overlap or less overlap between frames).

In a manner similar to that of comparative systems as discussed above with respect to FIG. 1B, some aspects of embodiments of the present disclosure relate to using the received far-end signal x(t) to estimate an acoustic echo. In more detail, in some embodiments, the far-end signal x(t) may also be applied to a far-end signal feature extraction module 220. In some embodiments, the far-end signal feature extraction module 220 is substantially similar to the microphone signal feature extraction module 210 and includes an STFT module 222, an absolute value module 224, and a logarithmic operation module 226. The far-end signal feature extraction module 220 computes far-end signal features {tilde over (X)}_(k,f) (in the feature space or hyperspace) from the far-end signal x(t). According to some embodiments of the present disclosure, the far-end signal features {tilde over (X)}_(k,f) are supplied to an echo estimator 230, which is configured to compute estimated echo features {tilde over (V)}_(k,f) (in the feature space or hyperspace).

In some embodiments of the present disclosure, a near-end estimator 250 accepts the microphone signal features {tilde over (D)}_(k,f), the far-end signal features {tilde over (X)}_(k,f), and the estimated echo features {tilde over (V)}_(k,f) (or another output of the echo estimator 230) to compute estimated near-end speech features {tilde over (Q)}_(k,f). The estimated near-end speech features {tilde over (Q)}_(k,f) may then be supplied to feature inversion module (or signal synthesis module) 270, which may include an exponential operation module 272 (to invert the logarithmic operation applied to the input signals) and an inverse short time Fourier transform (iSTFT) module 274 to transform the estimated near-end speech features {tilde over (Q)}_(k,f) from the feature space or hyperspace to a time domain signal q(t), which is an estimate of the near-end speech or near-end signal s(t).

In various speech processing applications, using past and/or future frames of data can help in computing estimates characteristics of the current frame. In some of such speech processing applications, a fixed context window is used as the input to a fully-connected first layer of a deep neural network. In these comparative methods, the contextual information can be lost after this first layer as the information flows through deeper layers.

Accordingly, some aspects of embodiments of the present disclosure use the context features for both inputs and outputs of the neural network in order to keep the contextual information available throughout the neural network. According to some embodiments, the input features for a current frame includes the feature vector {tilde over (X)}_(k) of current frame k and feature vectors ({tilde over (X)}_(k-1), {tilde over (X)}_(k-2) . . . {tilde over (X)}_(k-6)) of six previous frames or causal frames (k−1, k−2, . . . , k−6). According to some embodiments of the present disclosure, causal windows (using only data from previous frames, as opposed to future frames) are chosen to prevent extra latency (e.g., when using causal windows of frames there is no need to wait for the arrival of future frames k+1, k+2, . . . before processing a current frame k). The seven frames with 50% overlap of the embodiment discussed above creates a receptive filed of 112 ms, which is generally long enough for processing the speech signal. To incorporate context awareness, some aspects of embodiments of the present disclosure relate to the use of unrolled deep gated recurrent unit (GRU) networks with seven time-steps (or frames) for both the echo estimation module and the near-end estimation module. However, embodiments of the present disclosure are not limited thereto and may be implemented with more than six prior frames of data or fewer than six prior frames of data.

FIG. 3 is a block diagram illustrating an input layer of an unrolled deep gated recurrent unit (GRU) network, according to one embodiment of the present disclosure. As shown in FIG. 3, each of seven GRUs of the input GRU layer 232 receives a corresponding input feature vector from the feature vectors corresponding to the current frame {tilde over (X)}_(k) and the six previous frames {tilde over (X)}_(k-1), {tilde over (X)}_(k-2), . . . {tilde over (X)}_(k-6). Each of the seven GRUs of the input GRU layer 232 computes a hidden feature vector or activation h. For example, the GRU of the input GRU layer 232 corresponding to the k−6th frame receives the corresponding feature vector {tilde over (X)}_(k-6) and computes activation h_(k-6) ¹. The activation from each prior GRU of the input GRU layer 232 is supplied to the next GRU in the sequence. For example, activation {tilde over (X)}_(k-6) from the GRU corresponding to the k−6th frame is supplied as an input to the GRU for the k−5th frame. Accordingly, the GRU for the k−5th frame computes its activation h_(k-5) ¹ from its corresponding input feature vector {tilde over (X)}_(k-5) and the activation h_(k-6) ¹ from the previous frame. This recurrent arrangement allows contextual information from prior frames to be used when processing a current frame.

According to some embodiments of the present disclosure, each GRU computes its output activation in accordance with: h _(k) =z _(k) ⊙h _(k-1)+(1−z _(k))⊙ĥ _(k) where ⊙ is an element-wise multiplication, and the update gates z_(k) are: z _(k)=σ(W _(z) {tilde over (X)} _(k) +U _(z) h _(k-1)) where σ is a sigmoid function. The candidate hidden state ĥ_(k) is computed by ĥ _(k)=elu(W{tilde over (X)} _(k) +U(r _(k) ⊙h _(k-1))) where elu is exponential linear unit function, and reset gates r_(k) are computed by r _(k)=σ(W _(r) {tilde over (X)} _(k) +U _(r) h _(k-1)) where U, W, U_(r), W_(r), U_(z), and W_(z) are the internal weight matrices of the GRUs. In some embodiments, each of the GRUs in a given layer (e.g., each of the GRUs of the input GRU layer 232) uses the same set of weights (hence the “recurrent” nature of the neural network). In some embodiments, the values of the internal weight matrices are learned through a training process, described in more detail below.

FIG. 4 is a block diagram illustrating an architecture of an acoustic echo cancellation (AEC) neural network, according to one embodiment of the present disclosure. An acoustic echo cancellation (AEC) neural network 228 includes a deep neural network model that includes two stacked GRU networks. The first stack corresponds to the echo estimator 230 shown in FIG. 2 and takes context-aware frames {tilde over (X)}_(k) as inputs to each GRU of the input GRU layer 232 and estimates echo features {tilde over (V)}_(k) using a fully connected (FC) output layer 236 with linear activation. In the embodiment shown in FIG. 4, the first stack includes one hidden GRU layer 234. In more detail, the first GRU layer or input GRU layer 232 computes first activations h_(k-6) ¹, h_(k-5) ¹, . . . , h_(k) ¹, as discussed above with respect to FIG. 3. The first activations are supplied to the second GRU layer 234 to compute second activations h_(k-6) ², h_(k-5) ², . . . , h_(k) ². In the embodiment shown in FIG. 4, the second GRU layer 234 is also the last GRU layer of the stack. As such, the activations of the last GRU layer of the stack (here, the second activations h_(k-6) ², h_(k-5) ², . . . , h_(k) ²) are supplied to the fully connected output layer 236 to compute the estimated echo features {tilde over (V)}_(k-6), {tilde over (V)}_(k-5), . . . {tilde over (V)}_(k). However, embodiments of the present disclosure are not limited thereto and may include more than one hidden layer in the neural network of the echo estimator 230. For example, one or more additional GRU layers may be inserted between the input GRU layer 232 and the last GRU layer 234, where each i-th additional layer would compute its own activations h_(k-6) ^(i), h_(k-5) ^(i), . . . , h_(k) ^(i) based on the activations of the previous layer h_(k-6) ^(i-1), h_(k-5) ^(i-1), . . . , h_(k) ^(i-1).

In the embodiment shown in FIG. 2, the outputs h_(k-6) ², h_(k-5) ², . . . , h_(k) ² of the last GRU layer 234 from the first stack (the echo estimator 230) are supplied to a concatenation layer 251 to be concatenated with the context-aware frames corresponding to the features of the far-end signal {tilde over (X)}_(k) and the features {tilde over (D)}_(k) of the microphone signal d(t) to create concatenated inputs c to the GRU layers of the second stack (the near-end estimator 250). In the embodiment shown in FIG. 4, the second stack (the near-end estimator 250) includes three GRU layers: a third GRU layer 253 (which receives the concatenated inputs c), a fourth GRU layer 254, and a fifth GRU layer 255, which compute corresponding third activations h_(k-6) ³, h_(k-5) ³, . . . , h_(k) ³, fourth activations h_(k-6) ⁴, h_(k-5) ⁴, . . . , h_(k) ⁴, and fifth activations h_(k-6) ⁵, h_(k-5) ⁵, . . . , h_(k) ³. The activations from the last of the hidden GRU layers (in the embodiment shown in FIG. 4, fifth activations h_(k-6) ⁵, h_(k-5) ⁵, . . . , h_(k) ⁵ from the fifth GRU layer 255) are supplied to a fully connected output layer 259 to estimate the context-aware frames {tilde over (Q)}_(k) of features of estimated near-end speech. As shown in FIG. 4, the network may also compute the estimated-near-end speech of the prior six frames {tilde over (Q)}_(k-6), {tilde over (Q)}_(k-5), . . . {tilde over (Q)}_(k-1), which may be used for computing a loss during training in accordance with a loss function, as discussed in more detail below.

As noted above, in the embodiment shown in FIG. 4, data from the current frame k and the six previous frames k−6, k−5, . . . , k−1 are used to compute the features of the estimated near-end speech {tilde over (Q)}_(k). Also as discussed above, the embodiment shown in FIG. 4 makes use of a 257-point feature vector. Accordingly, altogether, the seven frames of the feature vectors {tilde over (X)}_(k) of the far-end signal x(t) have dimensions 7×257. In the embodiment shown in FIG. 4, the activations or outputs h of each of the GRU units is a vector of length 1000, although embodiments of the present disclosure are not limited thereto. Because the embodiment of FIG. 4 uses seven frames of data, the outputs or activations of each of the GRU layers 232 and 234 has dimensions of 7×1000. To match the shape of the input feature vectors, each of the output estimated echo features {tilde over (V)} has a length of 257, and therefore the output of the fully connected output layer 236 of the echo estimator has dimensions 7×257.

In the embodiment shown in FIG. 4, the microphone signal d(t) is supplied to substantially the same feature extractor as that used with the far-end signal x(t), and therefore the microphone signal features {tilde over (D)} for each frame are also feature vectors of length 257 values for each frame, such that the seven frames of microphone signal features {tilde over (D)} result in overall dimensions of 7×257. Each of the concatenation units of the concatenation layer 251 concatenates the output (having length 1000) of the last GRU layer 234 (e.g., output activations h²) of the first stack, the far-end signal features {tilde over (X)} (having length of 257 values), and the microphone signal features {tilde over (D)} (having length 257) of a corresponding frame, such that the output c of each of the concatenation units has a length of 1000+257+257=1514, and the output of the concatenation layer 251 has dimensions of 7×1514.

In the embodiment shown in FIG. 4, in a manner similar to that of the echo estimator 230, each GRU layer 253, 254, and 255 of the near-end estimator 250 produces an output activation h of dimensions 7×1000 (e.g., each GRU unit map compute an activation or activation vector having a length of 1000 values), and the final fully connected layer 259 produces feature vectors {tilde over (Q)} (having length 257) of estimated near-end speech of the current frame and the previous six frames, where the dimensions of the output of the fully connected layer 259 are 7×257.

FIG. 5A is a flowchart depicting a method for training a recurrent neural network of an acoustic echo cancellation module, according to one embodiment of the present disclosure. Generally, a supervised learning process of training a neural network involves computing various parameters of the neural network, such as weights and biases of neurons of the network, by adjusting the parameters to minimize a loss function between the output of the neural network and the “ground truth” data of a set of labeled training data for given set of inputs in the training data. The training may be performed by a computer system (including a processor and memory, and which may be accelerated using a vector processor such as a graphics processing unit) having access to training data, where the training data may be divided into a training set, a test set, and, in some embodiments, a validation set. Generally, the training set is used to learn the parameters of the network, the validation set is used to determine various hyper-parameters of the network (e.g., architectural parameters such as the number of neurons in each unit and the number of layers in the neural network), and the test set is used to evaluate the overall performance of the trained system.

In the particular domain of acoustic echo cancellation described here, the training data may include: far-end signals x(t); near-end signals s(t); and echo signals y(t). In some embodiments of the present disclosure, a method 500 for training a recurrent neural network of an acoustic echo cancellation module includes at 510, the computer system generates training data in a manner similar to that described in H. Zhang and D. Wang, “Deep Learning for Acoustic Echo Cancellation in Noisy and Double-Talk Scenarios,” in Proc. Annual Conference of the International Speech Communication Association, 2018, pp. 3239-3243. In some embodiments, the TIMIT dataset is used to generate the training data (see, e.g., F. Lamel, R. H. Kassel, and S. Seneff, “Speech database development: Design and analysis of the acoustic-phonetic corpus,” in Speech Input/Output Assessment and Speech Databases, 1989).

FIG. 5B is a flowchart of a method for generating training data, according to one embodiment of the present disclosure. In more detail, in some embodiments, the input data set may include recorded speech from a plurality of different human speakers or voices. A method 510 for generating training data includes at 511, a plurality of pairs of human speakers may be chosen (e.g., randomly) to be used as the near-end and far-end speakers. At 513, from each pair, three utterances of the far-end speaker are randomly chosen and concatenated to create a ground truth far-end signal x(t). At 515, for each pair, the far-end signal x(t) is transformed to generate a corresponding echo signal y(t) by simulating the effect of the acoustic path. FIG. 5C, discussed in more detail below, describes a method for generating a simulated ground truth echo signal y(t) according to one embodiment of the present disclosure.

At 517, each utterance of a near-end speaker of the pair is padded or extended to the same length as that of its corresponding far-end signal x(t) (e.g., for each concatenated far-end signal generated in accordance with the paired far-end human speaker) by filling zeroes before and after the utterance to have the same size as the far-end signal to generate ground truth near-end signals s(t). (Embodiments of the present disclosure are not limited thereto, and, in some embodiments, noise is added to the entire padded signal.) In some embodiments, more than one far-end signal x(t) and near-end signal s(t) pair is selected for each near-end far-end pair.

At 519, the computer system mixes (e.g., adds) the ground truth echo signals y(t) and the ground truth near-end signals s(t) computed for each pair to generate a corresponding training microphone signal d(t). For training mixtures, in some embodiments, the computer system generates the training microphone signals d(t) at 519 at signal to echo ratio (SER) level randomly chosen from {−6, −3, 0, 3, 6}dB by mixing the near-end speech signal and echo signal. The SER level is calculated on the double-talk period as:

${{SER}({dB})} = {10\log_{10}\frac{E\left\{ {s^{2}(t)} \right\}}{E\left\{ {y^{2}(t)} \right\}}}$

FIG. 5C is a flowchart of a method 515 for generating a simulated or ground truth echo signal y(t) according to one embodiment of the present disclosure in a manner similar to that described in S. Malik and G. Enzner, “State-space frequency-domain adaptive filtering for nonlinear acoustic echo cancellation,” IEEE Transactions on audio, speech, and language processing, vol. 20, no. 7, pp. 2065-2079, 2012. For the nonlinear model of acoustic path, at 515-1, the computer system applies hard clipping to simulate the power amplifier of loudspeaker (in one embodiment, x_(max) is set to 80% of the maximum volume of input signal):

${x_{clip}(t)} = \left\{ \begin{matrix} {- x_{\max}} & {{{if}\mspace{14mu}{x(t)}} < {- x_{\max}}} \\ {x(t)} & {{{if}\mspace{9mu}{{x(t)}}} \leq x_{\max}} \\ x_{m\;{ax}} & {{{if}\mspace{20mu}{x(t)}} > x_{\max}} \end{matrix} \right.$

At 515-3, to simulate the loudspeaker distortion, the computer system applies the a sigmoidal function such as:

${x_{nl}(t)} = {4\left( {\frac{2}{1 + {\exp\left( {{- a} \cdot {b(t)}} \right)}} - 1} \right)}$ where b(t)=1.5x_(clip)(t)−0.3x_(clip)(t)² and α=4 if b(t)>0 and α=0.5 otherwise.

According to one embodiment, at 515-5, a room impulse response (RIR) g(t) is randomly chosen from a set of RIRs, where the length of each of the RIRs is 512, the simulation room size is 4 meters×4 meters×3 meters, and a simulated microphone is fixed at the location of [2 2 1.5] meters (at the center of the room). A simulated loudspeaker is placed at seven random places with 1.5 m distance from the microphone. In some embodiments of the present disclosure, a plurality of different RIRs are also generated with different room sizes and different placements of the simulated microphone and/or simulated speaker.

In some embodiments, the RIRs are generated using an image method (see, e.g., J. B. Allen, D. A. Berkley, “Image method for efficiently simulating small-room acoustics,” The Journal of Acoustic Society of America, vol. 65, no. 4, pp. 943-950, 1979.) at reverberation time (T₆₀) of 200 ms. From the generated RIRs, in some embodiments some of the RIRs are used to generate the training data (e.g., may be randomly selected) while others are reserved to generate test data.

At 515-7, the output of sigmoidal function is convolved with the randomly chosen room impulse response (RIR) g(t) in order to simulate the acoustic transmission of the distorted (nonlinear) far-end signal x_(nl)(t) played through the loudspeaker in the room: y _(nl)(t)=x _(nl)(t)*g(t) where * indicates a convolution operation.

In some embodiments, a linear acoustic path y_(lin)(t) is simulated by only convolving the original far-end signal x(t) with the RIR g(t) to generate the echo signal, where nonlinearities such as clipping and loudspeaker distortion are not applied for this model: y _(lin)(t)=x(t)*g(t)

Referring back to FIG. 5A, at 520 the computer system computes feature vectors (ground truth near-end features {tilde over (S)}, ground truth far-end features {tilde over (X)}, ground truth echo features {tilde over (Y)}, and microphone features {tilde over (D)}) from respective parts of the training data (ground truth near-end signal s(t), ground truth far-end signal x(t), ground truth echo signal y(t), and the microphone signal d(t)) using the feature extractors as described above.

At 530, the computer system trains the neural network of the AEC 228 in accordance with the training data. In more detail, as discussed above, each of the GRUs computes its corresponding activation h from its inputs based on internal weight matrices U, W, U_(r), W_(r), U_(z), and W_(z). In addition, each of the fully connected units includes a plurality of internal weights W and biases b (e.g., applying an affine function of the form Wx+b) for mapping the inputs to the fully connected units to the outputs in feature space (e.g., STFT space).

Training the neural network involves learning the internal weights of the GRUs and the fully connected units such that the output feature vectors (estimated near-end features {tilde over (Q)} and estimated echo features {tilde over (V)}) are close to the ground truth feature vectors (ground truth near-end features {tilde over (S)} and ground truth echo features {tilde over (Y)}). The difference between the output feature vectors {tilde over (Q)} and {tilde over (V)} and the ground truth feature vectors {tilde over (S)} and {tilde over (Y)} may be measured using a loss function, representing how well the neural network, as configured with the current set of internal weights, approximates the underlying data.

In one embodiment, a mean absolute error (MAE) loss function is used for training the neural network. A mean absolute error is calculated between a ground-truth source (near-end signal s(t)) and a network estimated output (estimated near-end signal q(t)) in the feature domain (e.g., the STFT domain, as discussed above). Some embodiments use a weighted loss function that accounts for both the near-end signal s(t) and the echo path signal y(t) to compute the network weights. Accordingly, in one embodiment, the loss for a given frame k is computed based on the current frame and the previous six frames in accordance with:

${loss}_{k} = {{\beta{\sum\limits_{n = 0}^{m}{{{\overset{˜}{S}}_{k - n} - {\overset{˜}{Q}}_{k - n}}}_{1}}} + {\left( {1 - \beta} \right){\sum\limits_{n = 0}^{m}{{{\overset{˜}{Y}}_{k - n} - {\overset{˜}{V}}_{k - n}}}_{1}}}}$ where β is the weighting factor between the loss associated with the near-end signal and the loss associated with the echo signal, {tilde over (S)}_(i) corresponds to the ground truth near-end features for an i-th frame, {tilde over (Q)}_(i) corresponds to the estimated near-end features for the i-th frame, {tilde over (Y)}_(i) corresponds to the ground truth echo features for the i-th frame, and {tilde over (V)}_(i) corresponds to the estimated echo features for the i-th frame. In embodiments where m previous frames of data are used for context (e.g., a causal window of length m frames), the summations run from n=0 to m. For the sake of convenience, in the embodiments described in detail herein, m=6.

In some embodiments of the present disclosure, the weights are computed using gradient descent and backpropagation. In particular, the weights are iteratively adjusted based on the differences between the current output of the neural network and the ground truth. In some embodiments of the present disclosure, the models are trained using AMSGrad optimization (see, e.g., J. Reddi, S. Kale, and S. Kumar, “On the convergence of Adam and beyond,” in International Conference on Learning Representations (ICLR), 2018), and in particular the Adam variant (see, e.g., D. P. Kingma and J. L. Ba, “Adam: a method for stochastic optimization,” in International Conference on Learning Representations (ICLR), 2015) by setting β₁=0.9, β₂=0.999, and ∈=10⁻³ for 100 epochs, with a batch size of 100. In some embodiments, the weights of all layers are initialized with the Xavier method (see, e.g., X. Glorot, and Y. Bengio, “Understanding the difficulty of training deep feedforward neural networks,” in Proc. International Conference on Artificial Intelligence and Statistics, 2010, pp. 249-256) and with the biases initialized to zero. In some embodiments, L2 regularization for all the weights with a regularization constant of 0.000001 is used to prevent overfitting.

After training the weights of the neural network, the trained network may be tested using the test set of the training data to verify the accuracy of the network. As noted above, the test set may be formed using utterances from speakers who were not used in the training set and/or use RIRs and/or other distortions that were not present in the training set. Accordingly, the test set may be used to evaluate whether that the training process has trained a neural network to perform a generalized function for acoustic echo cancellation, rather than overfitting to the particular characteristics of the training data (e.g., removing acoustic echoes characteristic of the particular human speakers or RIRs of the training data).

After training the neural network and determining that the performance of the trained network is sufficient (e.g., based on the test set), the weights may be saved and used to configure a neural network running on an end-user device such as a smartphone or a tablet computer. In various embodiments of the present disclosure, the neural network of the acoustic echo cancellation module is implemented on at least one processor 1620 of the near-end device 10 (see, e.g., FIG. 11), where the processor may be: a general purpose central processing unit; a graphical processing unit (GPU); a field programmable gate array (FPGA); an neural processing unit (NPU) or neural network processor (NNP) (e.g., a processor having an architecture tailored to perform inference using a neural network); or a neuromorphic processor. For example, the parameters of the neural network (e.g., weights and biases) and the neural network architecture may be stored in non-transitory memory connected to the processor, where the processor performs inference using the network by loading the parameters and network architecture from memory. As another example, in the case of an FPGA, the FPGA may be configured in a non-transitory manner with the network architecture and the weights using a bitfile. Because the training process may be considered to be complete or stable, in some embodiments of the present disclosure, the end-user device may only operate the neural network in inference mode to compute the current estimated near-end features {tilde over (Q)} or estimated near-end signal q(t).

FIG. 6A is a block diagram illustrating an architecture of an acoustic echo cancellation (AEC) neural network operating in inference mode, according to one embodiment of the present disclosure. As shown in FIG. 6A, during inference, it is not necessary to compute the estimated near-end features for the prior frames {tilde over (Q)}_(k-6), {tilde over (Q)}_(k-5), . . . {tilde over (Q)}_(k-1), e.g., because, when the AEC neural network is used in the field, the near-end device 10 need only compute and transmit the current estimated frame to the far-end (e.g., the prior estimated near-end frames have already been transmitted to the far-end). (During training, the prior frames {tilde over (Q)}_(k-6), {tilde over (Q)}_(k-5), . . . {tilde over (Q)}_(k-1) were useful for providing additional information to the loss function.) Likewise, as shown in FIG. 6A, the fully connected output layer 236 of the first stack of the echo estimator 230 may be omitted from the inference mode neural network because the second stack of the near-end estimator 250′ in inference mode depends only on the output h_(k-6) ², h_(k-5) ², . . . , h_(k) ² of the last GRU layer 234 of the first stack of the echo estimator 230′.

FIG. 6B is a flowchart depicting a method 600 for computing an estimated near-end signal from a received far-end signal and a microphone signal using a neural network in inference mode, according to one embodiment of the present disclosure. As shown in FIG. 6B, at 602, the acoustic echo cancellation system 200 receives the far-end signal x(t) and, at 604, the far-end signal feature extraction module 220 extracts far-end features {tilde over (X)} from the far-end signal x(t). In addition, at 606 the first stack of the neural network, corresponding to the echo estimator 230′ computes an echo estimator output from the far-end features {tilde over (X)}. (In some embodiments, the echo estimator output corresponds to the output of the h_(k-6) ², h_(k-5) ², . . . , h_(k) ² of the last GRU layer 234 of the first stack of the echo estimator 230′.)

Similarly, at 612, the acoustic echo cancellation system 200 receives the microphone signal d(t) and, at 614, the microphone signal feature extraction module 210 extracts microphone signal features {tilde over (D)} from the microphone signal d(t).

At 620, the second stack of the neural network, corresponding to the near-end estimator 250′, computes estimated near-end features {tilde over (Q)} from the far-end features {tilde over (X)}, the echo estimator features (e.g., h), and the microphone features {tilde over (D)}. As shown in FIG. 6A, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, features from the context of prior frames (e.g., the six prior frames k−6, k−5, . . . , k−1) in addition to the current frame k are also supplied to the near-end estimator 250′. In particular, as shown in FIG. 6A, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, far-end features {tilde over (X)}, echo estimator outputs h, and microphone features {tilde over (D)} from the current frame k and the six prior frames k−6, k−5, . . . , k−1 (e.g., far-end features {tilde over (X)}_(k-6), {tilde over (X)}_(k-5), . . . , {tilde over (X)}_(k), echo estimator outputs h_(k-6) ², h_(k-5) ², . . . , h_(k) ², and microphone features {tilde over (D)}_(k-6), {tilde over (D)}_(k-5), . . . , {tilde over (D)}_(k)).

At 622, feature inversion module 270 of the acoustic echo cancellation system 200 computes an estimated near-end signal q(t) for the current frame from the estimated near-end features {tilde over (Q)} of the current frame. As noted above, the features (e.g., the far-end signal features {tilde over (X)}, the microphone features {tilde over (D)}, and the estimated near-end features {tilde over (Q)} may be in a feature space or hyperspace such as STFT space (e.g., spectral features or spectral domain). Accordingly, in some embodiments, the feature inversion module 270 transforms the estimated spectral features {tilde over (Q)} from the feature space to a time domain signal q(t) suitable for playback on a speaker at a far-end device. As shown in FIG. 2, the phase ∠D_(k,f) of the microphone signal d(t) may also be used by the inverse short-time Fourier transform (iSTFT) module 274 when computing the estimated near-end signal q(t).

FIG. 7 is a block diagram illustrating a hybrid system including the deep multi-task acoustic echo cancellation system, according to one embodiment of the present disclosure. In these embodiments, the multitask GRU network has the same architecture as described above. However, for the training of this network and during inference, the output of a frequency domain normalized least mean square (NLMS) filter 700 is used instead of the microphone signal d(t). The far-end signal x(t) remains as another input to the acoustic echo cancellation system 200 (including the multitask GRU network), and the output targets are still the near-end features {tilde over (S)} (as estimated as {tilde over (Q)}) and the echo features {tilde over (Y)} (as estimated as {tilde over (V)}).

Experimental Results Using Neural Network Techniques

To evaluate the performance of an acoustic echo cancellation system 200 as described above, experiments were performed using training data generated from the TIMIT dataset (see, e.g., F. Lamel, R. H. Kassel, and S. Seneff, “Speech database development: Design and analysis of the acoustic-phonetic corpus,” in Speech Input/Output Assessment and Speech Databases, 1989). In some embodiments of the present disclosure, seven utterances of near-end speakers were used to generate 3,500 training mixtures where each near-end signal was mixed with five different far-end signals. From the remaining 430 speakers, 100 pairs of speakers were randomly chosen as the far-end and near-end speakers. To generate 300 testing mixtures, the same procedure as described above, but with only three utterances of near-end speakers, where each near-end signal was mixed with one far-end signal. Therefore, the testing mixtures are from human speakers that were not part of the training set.

Perceptual Evaluation of Speech Quality (PESQ) scores of unprocessed test mixtures for linear and nonlinear models (no echo cancellation) are shown in Table 1. The unprocessed PESQ scores are calculated by comparing the microphone signal against near-end signal during the double-talk period.

TABLE 1 PESQ scores for unprocessed test mixtures in linear and nonlinear models of acoustic path Acoustic Testing SER (dB) Path Model 0 3.5 7 Linear 1.87 2.11 2.34 Nonlinear 1.78 2.03 2.26

In some instances, echo return loss enhancement (ERLE) was used to evaluate the echo reduction that is achieved by the acoustic echo cancellation system 200 according to embodiments of the present disclosure during the single-talk situations where only the echo is present, where ERLE is defined as:

${{ERLE}({dB})} = {10\log_{10}\frac{E\left\{ {d^{2}(t)} \right\}}{E\left\{ {q^{2}(t)} \right\}}}$ where E is the statistical expectation operation which is realized by averaging.

To evaluate the performance of the system during the double-talk periods, we used perceptual evaluation of speech quality (PESQ). In some embodiments, PESQ is calculated by comparing the estimated near-end speech q(t) against the ground-truth near-end speech s(t) during the double-talk only periods. A PESQ score ranges from −0.5 to 4.5 and a higher score indicates better quality.

In the following discussion, a frequency domain normalized least mean square (FDNLMS) (see, e.g., C. Faller and J. Chen, “Suppressing acoustic echo in a spectral envelope space,” IEEE Transactions on Acoustic, Speech and Signal Processing, vol. 13, no. 5, pp. 1048-1062, 2005) is used as a comparative example. A double-talk detector (DTD) is used based on the energy of microphone signal d(t) and far-end signal x(t). In some instances, a post-processing algorithm is further based on the method presented in R. Martin and S. Gustafsson, “The echo shaping approach to acoustic echo control”, Speech Communication, vol. 20, no. 3-4, pp. 181-190, 1996. Embodiments of the present disclosure are also compared against the bidirectional long short-term memory (BLSTM) method described in H. Zhang and D. Wang, “Deep Learning for Acoustic Echo Cancellation in Noisy and Double-Talk Scenarios,” in Proc. Annual Conference of the International Speech Communication Association, 2018, pp. 3239-3243.

Embodiments of the present disclosure are compared against comparative methods using a linear model of the acoustic path (e.g., linear acoustic echoes). Table 2 shows the average ERLE values and PESQ gains for the conventional NLMS filter, BLSTM, and a context-aware multitask GRU according to embodiments of the present disclosure (denoted as “CA Multitask GRU”). The PESQ gain is calculated as the difference of PESQ value of each method with respect to its unprocessed PESQ value. Table 2 also shows the results for context-aware single-task GRU (denoted as “CA Single-task GRU”) according to embodiments of the present disclosure that only uses the second stack of GRU layers with {tilde over (D)}_(k) and {tilde over (X)}_(k) as the inputs, where the loss function is calculated by only penalizing the network outputs against ground-truth feature vector {tilde over (S)} of near-end speech s(t). The results show that multitask GRU according to some embodiments of the present disclosure outperforms single-task GRU according to some embodiments of the present disclosure in terms of both PESQ and ERLE. It also shows that embodiments of the present disclosure outperform both conventional NLMS+Post-processing and BLSTM methods in all conditions.

TABLE 2 ERLE and PESQ scores in a linear model of acoustic path Testing SER (dB) Method 0 3.5 7 ERLE (dB) NLMS + Post-processing 29.38 25.88 21.97 BLSTM 51.61 50.04 47.42 CA Single-task GRU 62.88 61.81 60.11 CA Multitask GRU 64.66 64.16 62.26 PESQ gain NLMS + Post-processing 0.93 0.81 0.68 BLSTM 0.80 0.78 0.74 CA Single-task GRU 0.98 0.95 0.93 CA Multitask GRU 1.04 1.02 0.99

Embodiments of the present disclosure are also compared against comparative methods using a nonlinear model of the acoustic path (e.g., nonlinear acoustic echoes). In this set of experiments, the nonlinear ground truth echo signal y_(nl)(t) was used to generate the microphone signals d(t), therefore the model contains both power amplifier clipping and loudspeaker distortions (e.g., corresponding to 515-3 and 515-7 of FIG. 5C). The results of embodiments of the present disclosure are compared against comparative AES+RES methods including NLMS. In the nonlinear acoustic path, the performance was also compared against a deep neural network (DNN)-based residual echo suppression (RES) system that was described in C. M. Lee, J. W. Shin, and N. S. Kim, “DNN-based residual echo suppression,” in Proc. Annual Conference of the International Speech Communication Association, 2015, pp. 1775-1779. and denoted in the table as “AES+DNN.” The results in Table 3 show that embodiments of the present disclosure outperform the other two comparative methods in both PESQ and ERLE.

TABLE 3 ERLE and PESQ scores in nonlinear model of acoustic path Testing SER (dB) Method 0 3.5 7 ERLE (dB) NLMS + Post-processing 16.76 14.26 12.33 AES + DNN — 36.59 — CA Multitask GRU 61.79 60.52 59.47 PESQ gain NLMS + Post-processing  0.54 0.43  0.31 AES + DNN — 0.62 — CA Multitask GRU  0.84 0.83  0.81

Embodiments of the present disclosure achieve superior echo reduction without significant near-end distortion (e.g., the spectra corresponding to the estimated near-end signal and the actual near-end signal are very similar).

The performance of embodiments of the present disclosure was also evaluated in the presence of additive noise and a nonlinear model of the acoustic path. In these embodiments, when generating the training data, white noise at 10 dB SNR was added to the near-end signal s(t), with nonlinear acoustic path at 3.5 dB SER level. Embodiments of the present disclosure were then compared against a conventional NLMS+Post-processing system. As shown in Table 4 below, aspects of embodiments of the present disclosure outperform the comparative method by a large margin.

TABLE 4 ERLE and PESQ scores in nonlinear model of acoustic path (SER = 3.5 dB) and additive noise (SNR = 10 dB) ERLE (dB) NLMS + Post-processing 10.13 CA Multitask GRU 46.12 None 1.80 PESQ NLMS + Post-processing 2.01 CA Multitask GRU 2.50

In addition, the alternative hybrid embodiment discussed above was evaluated for unseen RIRs for different reverberation times and loudspeaker distances from the microphone. In this evaluation, the models were trained and tested using the same RIRs discussed above corresponding to a room size of 4 meters×4 meters×3 meters with reverberation time of 200 ms, and random loudspeaker distance of 1.5 meters from microphone and total length of 512 samples. During the testing of a hybrid system according to embodiments of the present disclosure, the loudspeaker distance was changed 15 cm. The results of frequency domain NLMS and a hybrid method of NLMS and multitask GRU according to embodiments of the present disclosure that was trained with the above RIRs are shown in Table 5. The multitask GRU was further fine-tuned with the RIRs that were generated in multiple room sizes (small, medium, and large), various reverberation times (from 250 ms to 900 ms), and loudspeaker distance of 15 cm. The fine-tuned results are also shown in Table 5, below. These results suggest that the hybrid method according to some embodiments of the present disclosure can perform better if the model is fine-tuned with the impulse response of the target device (e.g., target end-user near-end device).

TABLE 5 ERLE and PESQ scores of hybrid method ERLE (dB) NLMS 14.70 Hybrid Multitask GRU 37.68 Hybrid Multitask GRU (Fine-tuned) 41.17 None 2.06 PESQ NLMS 2.70 Hybrid Multitask GRU 3.23 Hybrid Multitask GRU (Fine-tuned) 3.37

Additional Embodiments

Some embodiments of the present disclosure are directed to different architectures for the neural network of the acoustic echo cancellation system 200. FIG. 8A is a schematic block diagram of an acoustic echo cancellation system in which a deep neural network (denoted as “EchoDNN”) is used to cancel echoes from a microphone signal d(t), according to one embodiment of the present disclosure. According to some embodiments of the present disclosure, the EchoDNN 900 uses only multiple fully connected (FC) layers.

FIG. 8B is an architectural block diagram of a deep neural network configured to cancel echoes from a microphone signal, according to one embodiment of the present disclosure. As shown in FIG. 8B, the extracted features {tilde over (D)} and {tilde over (X)} from both microphone d(t) and far-end x(t) signals are fed into an input layer 902 of the fully connected deep network (EchoDNN) 900. In the embodiment shown in FIG. 8B, the network includes three hidden layers 904, 906, and 908, and an output layer 910. Each of these layers may be a fully connected layer implementing, for example, an affine transformation. In some embodiments, each of three fully connected hidden layers 904, 906, and 908 include 4,000 units (e.g., neurons). In some embodiments, the output layer 910 is also a fully connected layer with 257 units (equal to the number of frequency bins of the extracted Log-Mag features). For the Log-Mel-Mag features, in some embodiments, the output layer used 80 units. In some embodiments, an exponential linear unit (ELU) is used as the activation function for each unit.

The estimated features 912 of the near-end signal are obtained directly from the output of the fully connected deep network (EchoDNN) 900. These features are converted back to the time-domain at 972 to synthesize the estimated near-end speech signal, e.g., using the feature inversion module 270 described above. In some embodiments, for both microphone d(t) and near-end signals s(t) sampled at the rate of 16 kHz, a frame size of 512 samples with 50% overlap was used. A 512-point short-time Fourier transform (STFT) was then applied to each frame of input signals resulted in 257 frequency bins. The final log-magnitude (Log-Mag) features were computed after calculating the logarithm operation on the magnitude values. In some embodiments of the present disclosure, the log-mel-magnitude (Log-Mel-Mag) was used as the final features 912 to reduce the dimensionality of the feature space and therefore reduce the complexity of the technique applied in these embodiments. In some embodiments, the features are compressed by using a 80-dimensional Mel-transformation matrix.

In order to use contextual information, in some embodiments, features for contextual frames of both input signals are also extracted and concatenated as the input features.

In various embodiments either log-magnitude (Log-Mag) features or Log-mag (or Log-Mel-Mag) of the near-end speech signal were used as the target labels during training.

In some embodiments, AMSGRAD is used as the optimizer during training. In some embodiments, mean absolute error (MAE) between the target labels and the output of the network was used as the loss function.

FIG. 9 is a block diagram illustrating an acoustic echo cancellation system using a generative network, according to one embodiment of the present disclosure. In the embodiment shown in FIG. 9, a generative model 990 such as WaveNet (see, e.g., Oord, Aaron van den, et al. “Wavenet: A generative model for raw audio.” arXiv preprint arXiv:1609.03499 (2016)) at the output of the AEC network (e.g., at the output of the near-end estimator 250) in place of the feature inversion module 270. Generative networks such as WaveNet can enhance the output of an acoustic echo cancellation system (AEC) 228 based on the information (e.g., mel-log-magnitude) that is provided by the AEC output or “Mel-spec.” This generative model 990 reconstructs the near-end speech signal q(t) at a much higher accuracy in comparison to the inverse STFT (iSTFT) 274.

Echo Removal Using a Contextual Attention Neural Network

Attention-based models provide one category of approaches to performing sequence-to-sequence learning tasks (e.g., a sequence of input audio frames to a sequence of output audio frames). Some aspects of embodiments of the present disclosure relate to contextual attention neural networks (e.g., multi-head attention networks) for acoustic echo cancellation (AEC). One example embodiment includes an encoder/decoder architecture with a multi-head contextual self-attention mechanism (or neural network) that maps the spectral features of microphone signals and far-end signals to estimate the near-end signal. Experiments using both simulated and real recorded room impulse responses (RIRs) show that contextual attention approaches according to some embodiments of the present disclosure achieve better performance in terms of the echo return loss enhancement (ERLE) score and the perceptual evaluation of speech quality (PESQ) score compared to benchmark echo removal techniques such as an encoder/decoder architecture lacking an attention-based component.

FIG. 10 is a block diagram illustrating an acoustic echo cancellation system including a contextual attention neural network, according to one embodiment of the present disclosure. The embodiment shown in FIG. 10 is generally similar to the embodiment shown in FIG. 2, but differs in that its echo canceller 1028 does not include a separate echo estimator (e.g., an echo estimator 230 similar to that shown in FIGS. 2 and 4). Instead, as described in more detail below, the echo canceler 1028 includes a contextual attention neural network configured to compute the estimated near end speech features {tilde over (Q)}_(k,f) from the microphone signal features {tilde over (D)}_(k,f) and the far-end signal features {tilde over (X)}_(k,f) without training a network to also estimate echo features {tilde over (V)}_(k,f). Like reference numerals in FIG. 10 correspond to like components of the embodiment shown in FIG. 2 and therefore the similar components of FIG. 10 will not be described in detail.

As shown in FIG. 10, like in the embodiment of FIG. 2, the microphone signal d(t)=s(t)+y(t) is supplied to the acoustic echo cancellation system 1000 to cancel or remove the echo signal y(t). The microphone signal d(t) is supplied to a microphone signal feature extraction module 1010, which includes a short time Fourier transform (STFT) module 1012, an absolute value module 1014 to remove the conjugate symmetric half of the output of the STFT module 1012, and a logarithmic operation module 1016 to compute the final logarithmic magnitude spectral feature vector or logarithmic spectral features or microphone signal features {tilde over (D)}_(k,f) in “feature space” or “hyperspace” or “logarithmic spectral space.” The STFT module 1012 may also compute a phase ∠D_(k,f) of the microphone signal d(t).

As shown in FIG. 10, like in the embodiment of FIG. 2, the far-end signal x(t) may also be applied to a far-end signal feature extraction module 1020. In some embodiments, the far-end signal feature extraction module 1020 is substantially similar to the microphone signal feature extraction module 1010 and includes an STFT module 1022, an absolute value module 1024, and a logarithmic operation module 1026. The far-end signal feature extraction module 1020 computes far-end signal features {tilde over (X)}_(k,f) (in the feature space or hyperspace) from the far-end signal x(t).

In the embodiment shown in FIG. 10, the microphone signal features {tilde over (D)}_(k,f) and the far-end signal features {tilde over (X)}_(k,f) are supplied to an echo canceller 1028 (or acoustic echo canceller), which includes a contextual attention neural network (which may be referred to as a contextual attention aware model).

FIG. 11A is a block diagram depicting a contextual attention neural network, according to one embodiment of the present disclosure. As noted above, in some embodiments of the present disclosure, a contextual attention neural network 1100 is a component of the echo canceller 1028 shown in FIG. 10 and is used in the process of computing near-end features {tilde over (Q)}_(k,f) based on 10, the microphone signal features {tilde over (D)}_(k,f) and the far-end signal features {tilde over (X)}_(k,f).

FIG. 11B is a more detailed block diagram depicting layers of gated recurrent units (GRUs) of a contextual attention neural network, according to one embodiment of the present disclosure.

As shown in FIG. 11A and FIG. 11B, the contextual attention neural network 1100 includes an encoder module 1110, a contextual attention module 1130, and a decoder module 1190, where each of the modules includes one or more layers of neurons (e.g., layers of gated recurrent units), as discussed in more detail below. The encoder module 1110 is configured to receive the microphone signal features {tilde over (D)}_(k,f) and the far-end signal features {tilde over (X)}_(k,f) as input. In the embodiments shown in FIG. 11A, each of these two inputs has a dimension of T×257, where T is the number of frames of input used in the causal window or contextual window. For example, in the particular example embodiments shown and described above with respect to FIGS. 3, 4, and 6A, the number of frames T was set to seven frames (T=7), where a current frame and the six previous causal frames were used to perform the echo removal. As noted above, embodiments of the present disclosure are not limited to the case where the window has a size of seven frames (T=7). As noted above, embodiments of the present disclosure are described herein where the features are represented as 257-point vectors. Accordingly, as shown in FIG. 11A and FIG. 11B, each frame of the microphone signal features {tilde over (D)}_(k,f) and the far-end signal features {tilde over (X)}_(k,f) includes 257×2=514 features, and therefore the input to the network 1100 has a total dimension of T×514.

As shown in FIG. 11A, the encoder module 1110 takes the logarithmic magnitudes of the microphone signal features {tilde over (D)}_(k,f)=[{tilde over (D)}_(t-T+1), . . . , {tilde over (D)}_(t)] and the far-end signal features {tilde over (X)}_(k,f)=[{tilde over (X)}_(t-T+1), . . . {tilde over (X)}_(t)] and maps them to a hyperspace or latent space to generate encoded features h, where the encoded features in the embodiment shown in FIG. 11A have dimensions of T×1000: h=Encoder({tilde over (D)} _(k,f) ,{tilde over (X)} _(k,f))

In more detail, in some embodiments of the present disclosure as shown in FIG. 11B, the encoder module 1110 includes a concatenation layer 1112 configured to concatenate the input microphone signal features {tilde over (D)}_(k,f) and far-end signal features {tilde over (X)}_(k,f) (each having dimensions T×257) to compute concatenated features of dimension T×514. The concatenated features are then supplied to a gated recurrent unit (GRU) layer 1114 substantially similar to the GRU layer shown in FIG. 3, where each unit of the GRU layer 1114 computes the encoded features h in accordance with: h _(t) =f([{tilde over (D)} _(t) ;{tilde over (X)} _(t)],h _(t-1)) where f represents the GRU function and h_(t) represents the output of the GRU for time t.

The contextual attention module 1130 takes the encoded features h to compute contextual attentions (or contextual attention features) c that identify important regions of the hyperspace, where the contextual attentions c in the embodiment shown in FIG. 11A have dimensions of T×1000: c=Attention(h)

In more detail, the contextual attention mechanism allows the model to learn the importance of each speech frame within the context. This is done by weighting the speech frames when constructing the representation of the frames in the hyperspace or latent space. Applying contextual attention in accordance with embodiments of the present disclosure is particularly effective in double-talk periods. As noted above, while the attention mechanism may be applied to both past frames and future frames, in the interest of avoiding or reducing latency, some embodiments of the present disclosure use only past frames.

FIG. 11C is a block diagram depicting a contextual attention module, according to one embodiment of the present disclosure. The contextual attention module 1130 shown in FIG. 11C includes two multi-head attention (MHA) layers (or modules) 1131 and 1134, where the first multi-head attention layer 1131 implements a contextual self-attention layer. In general, contextual self-attention layers, such as the multi-head self-attention layer 1131, are used to capture the similarity of any hyperspace representation of a speech frame with respect to all the past neighboring speech frames in an input sequence.

For the sake of convenience, the computations performed by a multi-head self-attention (MHSA) layer will be described in more detail below in the context of the multi-head self-attention layer 1131. In some embodiments, the multi-head attention (MHA) module 1134 is architecturally substantially the same as the MHSA module 1131 but uses different inputs than the MHSA module 1131, as shown in FIG. 11C, and as described in more detail below. Weight parameters α_(t,t′) capture the similarity between the hidden state representations supplied as inputs to the layers (e.g., h_(t) and h_(t′) of concatenated speech frame logarithmic magnitudes [{tilde over (D)}_(t); {tilde over (X)}_(t)] and [{tilde over (D)}_(t′); {tilde over (X)}_(t′)]) at time steps t and t′ respectively. In some embodiments, the contextual self-attention mechanism computes the weight parameters α_(t,t′) in accordance with:

$a_{t,t^{\prime}} = \frac{\exp\left( e_{t,t^{\prime}} \right)}{\sum\limits_{i}{\exp\left( {e_{t},i} \right)}}$ where e_(t,t′) are attention importance scores computed for each time step. In some embodiments of the present disclosure, the attention importance scores e_(t,t′) are calculated using an “additive” method: e _(t,t′) =W _(α) ^(T) tan h(W _(t) h _(t) +W _(t′) h _(t′)) In some embodiments of the present disclosure, the attention importance scores e_(t,t′) are calculated using a “general” method: e _(t,t′) =h _(t) ^(T) W _(α) h _(t′)

In some embodiments, the attention importance scores e_(t,t′) are normalized after they are computed to construct a probability distribution over the speech frames (e.g., by applying a softmax activation function to the scores).

In some embodiments, the output of the multi-head self-attention layer 1131 is an attention-focused hidden state representation of a frame at time-step t given by the weighted summation of the hidden state representation h_(t′) of all past frames at time steps t′ and their similarities α_(t,t′) to the hidden state representation h_(t) of the current frame. More informally, the output of the multi-head self-attention layer 1131 represents how much to attend to a speech frame at any time-step, based on the neighborhood context: Σ_(t′)α_(t,t′) h _(t′)

To improve the general method for calculating the importance scores e_(t,t′), in some embodiments, a temperature parameter is used in the softmax activation function and three learnable weights (W_(a), W_(b), W_(e)) are used instead of just one weight W_(a):

${softmax}\left( \frac{\left( {hW}_{b} \right)\left( {hW}_{a} \right)^{T}}{\sqrt{d_{h}}} \right){hW}_{c}$ (see, e.g., G. Hinton, O. Vinyals, and J. Dean, “Distilling the Knowledge in a Neural Network,” in NIPS Deep Learning and Representation Learning Workshop, 2015).

In some embodiments of the present disclosure, as shown in FIG. 11C, the multi-head self-attention layer is used with a multi-head attention layer 1134 (see, e.g., A. Vaswani et al., “Attention is All you Need,” in Advances in Neural Information Processing Systems, I. Guyon, U. V Luxburg, S. Bengio, H. Wallach, R. Fergus, S. Vishwanathan, and R. Garnett, Eds. Curran Associates, Inc., 2017, pp. 5998-6008), where both the multi-head self-attention layer 1131 and the multi-head attention layer 1134 have skip connections and layer normalizations.

In more detail, as shown in FIG. 11C, a contextual attention module 1130 (or contextual attention layer) according to one embodiment includes a multi-head self-attention module 1131, which receives three copies of the encoded features h. The output of the multi-head self-attention module 1131 is supplied to a first summation node 1132, which adds the output of the multi-head self-attention module 1131 to another copy of the encoded features h. The output of the first summation node 1132 is then supplied to a first layer normalization module 1133 (see, e.g., Ba, Jimmy Lei, Jamie Ryan Kiros, and Geoffrey E. Hinton. “Layer normalization.” arXiv preprint arXiv: 1607.06450 (2016)), and the layer normalized activations are supplied as an input to a multi-head attention module 1134 along with two additional copies of the encoded features h. In other words, the query for the multi-head attention layer 1134 is the output of the first layer normalization module 1133 and the output of the encoder (the encoded features h). The output of the multi-head attention module 1134 is supplied to a summation node 1135, which adds a copy of the output of the first layer normalization module 1133 before outputting the sum to a second layer normalization layer 1136. The output of the second layer normalization layer is then output from the contextual attention module 1130 as the contextual attentions (or contextual attention features) c.

In some embodiments of the present disclosure, the multi-head attention module 1134 and the second layer normalization module 1136 are omitted, such that the output of the first layer normalization module 1133 is taken as the output c of the contextual attention module 1130.

Referring back to FIGS. 11A and 11B, the estimated logarithmic magnitude of the near-end speech {tilde over (Q)}_(k,f) is then computed from the contextual attentions c using the decoder module 1190, where the dimensions of the near-end speech {tilde over (Q)}_(k,f) in the embodiments shown in FIG. 11A and FIG. 11B have dimensions of T×257 (e.g., the same dimensions as the input microphone and far-end signal features).

In more detail, in some embodiments, such as that shown in FIG. 11B, the decoder includes two GRU layers 1192 and 1194. According to one embodiment, the first GRU layer 1192 computes its outputs s (e.g., s_(k-T+1), s_(k-T+2), . . . s_(k)) based on the contextual attentions c (e.g., c_(k-T+1), c_(k-T+2), . . . , c_(k)) from the contextual attention module 1130 in accordance with: s _(t) =g ₁(c _(t) ,s _(t-1)) and the second GRU layer 1194 computes the spectral features of the estimated near end speech {tilde over (Q)}_(k,f) (e.g., {tilde over (Q)}_(k-T+1), {tilde over (Q)}_(k-T+2), . . . , {tilde over (Q)}_(k)) in accordance with: {tilde over (Q)} _(t) =g ₂(s _(t) ,{tilde over (Q)} _(t-1))

In some embodiments, an exponential linear unit (elu) activation function is used with the first GRU layer 1192 and a linear activation function is used with the second GRU layer 1194.

As shown in FIG. 10, the estimated features of the near-end speech near-end speech {tilde over (Q)}_(k,f) are converted back into estimated near-end speech signal q(t), which may be transmitted to the far-end communication device.

In some embodiments, the contextual attention neural network 1100 shown in FIGS. 11A, 11B, and 11C is trained using substantially the same training procedure as that described above with respect to FIGS. 5A, 5B, and 5C, and the resulting trained contextual attention neural network 1100 may be deployed into a communication device (e.g., an end-user device such as a smartphone or a tablet computer) in a manner similar to that described above.

FIG. 12 is a flowchart depicting a method 1200 for computing an estimated near-end signal from a received far-end signal and a microphone signal using a contextual attention neural network 1100 in inference mode, according to one embodiment of the present disclosure. As shown in FIG. 12, at 1202, the acoustic echo cancellation system 1000 receives the far-end signal x(t) and, at 1204, the far-end signal feature extraction module 1020 extracts far-end features {tilde over (X)} from the far-end signal x(t).

Similarly, at 1206, the acoustic echo cancellation system 1000 receives the microphone signal d(t) and, at 1208, the microphone signal feature extraction module 1010 extracts microphone signal features {tilde over (D)} from the microphone signal d(t).

At 1210, the contextual attention neural network 1100 (e.g., a component of the echo canceller 1028) computes estimated near-end features {tilde over (Q)} from the far-end features {tilde over (X)}, and the microphone features {tilde over (D)}. As shown, for example, in FIG. 11B, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, features from the context of prior T frames (e.g., the six prior frames k−6, k−5, . . . , k−1) in addition to the current frame k are also supplied to the contextual attention neural network 1100. In particular, as shown in FIG. 11B, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, far-end features {tilde over (X)} and microphone features {tilde over (D)} from the current frame k and the T−1 prior frames k−T+1, k−T+2, . . . , k−1 (e.g., far-end features {tilde over (X)}_(k-T+1), {tilde over (X)}_(k-T+2), . . . , {tilde over (X)}_(k) and microphone features {tilde over (D)}_(k-T+1), {tilde over (D)}_(k-T+2), . . . . , {tilde over (D)}_(k)).

At 1212, feature inversion module 1070 of the acoustic echo cancellation system 1000 computes an estimated near-end signal q(t) for the current frame from the estimated near-end features {tilde over (Q)} of the current frame. As noted above, the features (e.g., the far-end signal features {tilde over (X)}, the microphone features {tilde over (D)}, and the estimated near-end features {tilde over (Q)} may be in a feature space or hyperspace such as STFT space (e.g., spectral features or spectral domain). Accordingly, in some embodiments, the feature inversion module 1070 transforms the estimated spectral features {tilde over (Q)} from the feature space to a time domain signal q(t) suitable for playback on a speaker at the far-end communication device, where the feature inversion module 1070 may include an exponential module 1072 and an inverse short-time Fourier transform (iSTFT) module 1074. As shown in FIG. 10, the phase ∠D_(k,f) of the microphone signal d(t) may also be used by the inverse short-time Fourier transform (iSTFT) module 1074 when computing the estimated near-end signal q(t).

Experimental Results of Echo Cancellation Using Contextual Attention Neural Network Techniques

Embodiments of acoustic echo cancellation systems 1000 using a contextual attention neural network 1100 for removing acoustic echoes in accordance with embodiments of the present disclosure were trained and tested in a manner similar to that described above for testing the acoustic echo cancellation system 200.

To create a benchmark for comparison, the training data generated from the TIMIT data set was used to measure the performance of a frequency domain normalized least mean square (FDNLMS) model with double talk detection (DTD) based energy of the microphone and far-end signals (see, e.g., C. Faller and J. Chen, “Suppressing acoustic echo in a spectral envelope space,” IEEE Trans. Speech Audio Process., vol. 13, no. 5, pp. 1048-1061, 2005).

Some additional tests were performed on the effect of applying a post-filtering method to the FDNMLS audio (FDNMLS+Post-filtering), based on R. Martin and S. Gustafsson, “The echo shaping approach to acoustic echo control,” Speech Commun., vol. 20, no. 3-4, pp. 181-190, 1996. Further tests were performed using a deep neural network (DNN) instead of post-filtering, using a method as described in C. M. Lee, J. W. Shin, and N. S. Kim, “DNN-based residual echo suppression,” in INTERSPEECH, 2015, vol. 2015-January, pp. 1775-1779. where the parameters of the DNN were set based on the values given in the paper. (FDNLMS+DNN).

A fourth benchmark used substantially the same encoder/decoder GRU network architecture shown in FIG. 11A, but excluded the contextual attention module 1130 (e.g., the output of the encoder module 1110 was supplied directly as input to the decoder module 1190). (Encoder/Decoder GRU).

Embodiments of the preset disclosure corresponding to five different variations for computing attention were as described above. These include the use of general, additive, or scaled dot-product attention mechanisms, as well as the use of multi-head self-attention (MHSA) alone and the use of MHSA with multi-head attention (MHSA+MHA). Table 6 shows ERLE and PESQ scores in a linear model of acoustic paths using synthetic RIRs

TABLE 6 ERLE and PESQ scores in a linear model of acoustic path Testing SER (dB) Method 0 3.5 7 ERLE FDNLMS 14.20 13.10 11.72 (dB) FDNLMS + Post-filtering 29.38 25.88 21.97 FDNLMS + DNN 35.08 35.35 36.86 Encoder/Decoder GRU 63.82 63.26 61.17 Contextual Attention GRU General 62.45 62.97 61.73 Additive 63.30 63.51 62.01 Scaled Dot-Product 61.39 62.45 61.86 MHSA 63.44 62.78 61.53 MHSA + MHA 64.16 63.60 62.07 PESQ Unprocessed 1.87 2.11 2.34 FDNLMS 2.46 2.65 2.83 FDNLMS + Post-filtering 2.80 2.92 3.02 FDNLMS + DNN 2.56 2.76 2.94 Encoder/Decoder GRU 2.80 3.02 3.20 Contextual Attention GRU General 2.83 3.05 3.22 Additive 2.90 3.13 3.31 Scaled Dot-Product 2.91 3.13 3.32 MHSA 2.93 3.14 3.32 MHSA + MHA 2.94 3.16 3.34

As shown in Table 6, methods according to embodiments of the present disclosure using contextual attention generally outperform the benchmark echo removal systems, including FDNLMS, FDNLMS+Post-filtering, FDNLMS+DNN, and the bare encoder/decoder GRU without a contextual attention module. Among the variations of embodiments of the present disclosure, the use of MHSA and MHA together demonstrated the best ERLE and PESQ scores.

Embodiments of the present disclosure were also compared against the FDNLMS+Post-filtering and Encoder/Decoder GRU benchmarks in removing echoes due to non-linear acoustic paths, including simulated power amplifier clipping and loudspeaker distortions. Table 7 shows ERLE and PESQ scores in a non-linear model of acoustic path (SER=3.5 dB) using synthetic RIRs.

TABLE 7 ERLE and PESQ scores in a non-linear model of acoustic path ERLE FDNLMS + Post-filtering 14.26 (dB) Encoder/Decoder GRU 50.84 Contextual Attention GRU 57.06 PESQ Unprocessed 2.03 FDNLMS + Post-filtering 2.46 Encoder/Decoder GRU 2.84 Contextual Attention GRU 2.94

As shown in Table 7, above, a Contextual Attention GRU echo canceller in accordance with an embodiment of the present disclosure outperforms the FDNLMS+Post-filtering and Encoder/Decoder GRU (without a contextual attention module) on both ERLE and PESQ scores when experimenting with synthetic RIRs.

Table 8, below, further shows ERLE and PESQ scores of a Contextual Attention GRU echo canceller in comparison against the same FDNLMS+Post-filtering and Encoder/Decoder GRU (without a contextual attention module) baseline models using real (as opposed to synthetic) recorded RIRs from the Aachen impulse response database (see, e.g., M. Jeub, M. Schäfer, and P. Vary, “A binaural room impulse response database for the evaluation of dereverberation algorithms,” in International Conference on Digital Signal Processing, 2009).

TABLE 8 ERLE and PESQ scores in real recorded room impulse responses Testing SER (dB) Method 0 3.5 7 ERLE FDNLMS + Post-filtering 21.42 20.02 18.08 Encoder/Decoder GRU 55.12 58.80 60.23 Contextual Attention GRU 55.54 59.36 61.72 PESQ Unprocessed 1.86 2.10 2.33 FDNLMS + Post-filtering 2.71 2.88 3.03 Encoder/Decoder GRU 2.74 2.96 3.15 Contextual Attention GRU 2.88 3.10 3.28

As shown in Table 8, above, a Contextual Attention GRU echo canceller in accordance with an embodiment of the present disclosure outperforms the FDNLMS+Post-filtering and Encoder/Decoder GRU (without a contextual attention module) on both ERLE and PESQ scores when testing on real recorded RIRs.

Accordingly, some embodiments of the present disclosure are directed to the use of a contextual attention neural network as a component of an echo removal system. Such contextual attention neural networks provide increased performance (e.g., reduce the perceived echo in the transmitted near-end signal) over comparative echo removal techniques.

Echo Removal Using Frequency Domain Normalized Mean Square (FDNLMS) with Contextual Attention Neural Networks

Some aspects of embodiments of the present disclosure relate to the use of Frequency Domain Normalized Mean Square (FDNLMS) in conjunction with contextual attention neural networks similar to those described above.

In some comparative acoustic echo cancellation (AEC) techniques, the frequency domain normalized mean square (FDNLMS) is used to estimate the acoustic path. The received far-end signal is then modified by the acoustic path estimated using FDNLMS to compute an estimated echo signal, which is subtracted from the microphone signal.

The embodiments of the present disclosure described in this section relate to a hybrid approach that combines a classical echo estimation technique (FDNLMS is used as an example herein) with recurrent deep neural networks for AEC. In more detail, some aspects of embodiments of the present disclosure use the far-end signal and the microphone signal, along with the error signal between them, in the feature domain, to estimate the near-end signal using an encoder/decoder gated recurrent unit (GRU) networks. In some embodiments, the encoder/decoder GRU networks further include a contextual attention module between the encoder and decoder modules. The FDNLMS update rule is used to calculate the echo estimate. Using a hybrid approach combines the advantages of classical acoustic echo cancellation (e.g., dealing well with previously unseen environments) with the superior performance of deep learning methods. Experiments using both simulated and real recorded room impulse responses (RIRs) show that hybrid approaches (e.g., using FDNLMS with encoder/decoder GRU networks) in accordance with embodiments of the present disclosure can consistently achieve better performance in terms of the echo return loss enhancement (ERLE) and the perceptual evaluation of speech quality (PESQ) scores compared to various baseline classical and deep learning techniques.

FIG. 13 is a block diagram illustrating an acoustic echo cancellation system including a hybrid of a Frequency Domain Normalized Mean Square (FDNLMS) technique in conjunction with a contextual attention neural network, according to one embodiment of the present disclosure. The embodiment shown in FIG. 13 is generally similar to the embodiment shown in FIG. 2, but differs in that its echo canceller 1328 takes, as further input, error features {tilde over (E)}_(k,f) calculated in accordance with an FDNLMS technique. According to some embodiments of the present disclouse, the echo canceller 1328 includes a recurrent neural network. As described in more detail below, in some embodiments, the recurrent neural network of the echo canceller 1328 is a contextual attention neural network configured to compute the estimated near end speech features {tilde over (Q)}_(k,f) from the microphone signal features {tilde over (D)}_(k,f), the far-end signal features {tilde over (X)}_(k,f), and error signal features {tilde over (E)}_(k,f). Like reference numerals in FIG. 13 correspond to like components of the embodiments shown in FIG. 2 and FIG. 10 and therefore the similar components of FIG. 13 will not be described in detail.

As shown in FIG. 13, like in the embodiments of FIG. 2 and FIG. 10, the microphone signal d(t)=s(t)+y(t) is supplied to the acoustic echo cancellation system 1300 to cancel or remove the echo signal y(t). The microphone signal d(t) is supplied to a microphone signal feature extraction module 1310, which includes a short time Fourier transform (STFT) module 1312, an absolute value module 1314 to remove the conjugate symmetric half of the output of the STFT module 1312, and a logarithmic operation module 1316 to compute the final logarithmic magnitude spectral feature vector or logarithmic spectral features or microphone signal features {tilde over (D)}_(k,f) in “feature space” or “hyperspace” or “logarithmic spectral space.” The STFT module 1312 may also compute a phase ∠D_(k,f) of the microphone signal d(t).

As shown in FIG. 13, as in the embodiments of FIG. 2 and FIG. 10, the far-end signal x(t) may also be applied to a far-end signal feature extraction module 1320. In some embodiments, the far-end signal feature extraction module 1320 is substantially similar to the microphone signal feature extraction module 1310 and includes an STFT module 1322, an absolute value module 1324, and a logarithmic operation module 1326. The far-end signal feature extraction module 1320 computes far-end signal features {tilde over (X)}_(k,f) (in the feature space or hyperspace) from the far-end signal x(t).

FIG. 13 further depicts a classical echo estimator 1330 configured to compute error features {tilde over (E)}_(k,f). According to one embodiment of the present disclosure, the logarithmic spectral error features {tilde over (E)}_(k,f) are computed from the far-end signal x(t), the microphone signal d(t), and estimated acoustic path, which is estimated based on Frequency Domain Normalized Mean Square (FDNLMS) in the embodiment shown in FIG. 13. In more detail, absolute error features |E_(k,f)| may be computed in accordance with: |E _(k,f) |=|D _(k,f) |−G _(k,f) |X _(k,f)| where |D_(k,f)| are the absolute microphone features (e.g., the output of the absolute value module 1314 of the microphone signal feature extraction module 1310), |X_(k,f)| are the absolute far-end features (e.g., the output of the absolute value module 1324 of the far-end signal feature extraction module 1320), and G_(k,f) are parameters representing the current estimate of the acoustic path, as estimated using classical techniques. As shown in FIG. 13, a product node 1332 of the classical echo estimator 1330 computes the product G_(k,f)|X_(k,f)| of the acoustic path parameters G_(k,f) and the absolute far-end features |X_(k,f)|. This product represents an estimate of the echo signal y(t) and is depicted in FIG. 13 as absolute estimated echo features |{tilde over (Y)}_(k,f)|. When using FDNLMS, the acoustic path parameters G_(k,f) are computed and updated for each frame in accordance with:

$G_{{k + 1},f} = {G_{k,f} + {\frac{\mu}{P_{k,f}}{E_{k,f}}{X_{k,f}}}}$ where the step size μ is normalized by the average power P_(k,f) of the far-end signal x(t) and is obtained recursively by: P _(k,f)=(1−α)P _(k-1,f) +α|X _(k,f)|² where α is a forgetting factor between 0 and 1.

As shown in FIG. 13, a summation node 1334 subtracts the absolute estimated echo features |{tilde over (Y)}_(k,f)| from the absolute microphone features |D_(k,f)| to compute the absolute error features |E_(k,f)|, and a logarithmic operation module 1336 computes logarithmic error features {tilde over (E)}_(k,f) from the absolute error features |E_(k,f)|.

In the embodiment shown in FIG. 13, the microphone signal features {tilde over (D)}_(k,f) the far-end signal features {tilde over (X)}_(k,f), and the error features {tilde over (E)}_(k,f) are supplied to an echo canceller 1328. In some embodiments of the present disclosure, the echo canceller 1328 includes a recurrent neural network. As noted above, in some embodiments, the recurrent neural network is, or includes, a contextual attention aware model or a contextual attention neural network 1400, as shown, for example, in FIG. 14A and FIG. 14B. In some embodiments of the present disclosure, instead of using a contextual attention neural network 1400, the recurrent neural network is, or includes, includes an encoder/decoder GRU network without a contextual attention module (e.g., where the encoder supplies its output directly to the decoder).

FIG. 14A is a block diagram depicting a contextual attention neural network, where the contextual attention neural network is configured to accept far-end features, error features, and microphone features as input, according to one embodiment of the present disclosure. A contextual attention neural network 1400 may be used as a component of the echo canceller 1328 shown in FIG. 13 and is configured to accept far-end features {tilde over (X)}_(k,f), error features {tilde over (E)}_(k,f), and microphone features {tilde over (D)}_(k,f) as input. FIG. 14B is a more detailed block diagram depicting layers of a contextual attention neural network configured to accept far-end features, error features, and microphone features as input, according to one embodiment of the present disclosure. The contextual attention neural network 1400 is configured to accept far-end features {tilde over (X)}_(k,f), error features {tilde over (E)}_(k,f), and microphone features {tilde over (D)}_(k,f) as input. The contextual attention neural network 1400 shown in FIGS. 14A and 14B is substantially similar to that of the contextual attention neural network 1100 shown in FIGS. 11A and 11B, and therefore a description of the substantially components will not be repeated herein.

One difference from the contextual attention neural network 1100 shown in FIGS. 11A and 11B is that the contextual attention neural network 1400 takes as input the microphone signal features {tilde over (D)}_(k,f), the far-end signal features {tilde over (X)}_(k,f), and the error features {tilde over (E)}_(k,f) discussed above. In the examples shown in FIGS. 14A and 14B, each of the microphone signal features {tilde over (D)}_(k,f) the far-end signal features {tilde over (X)}_(k,f), and the error features {tilde over (E)}_(k,f) has dimensions T×257. Accordingly, when these features are concatenated (e.g., by a concatenation layer 1412 of an encoder module 1410 of the contextual attention neural network 1400), the resulting features have dimensions T×771.

In a manner similar to that of the embodiments shown in FIG. 11A and FIG. 11B, the contextual attention neural networks 1400 of the embodiments shown in FIGS. 14A and 14B include an encoder module 1410, a contextual attention module 1430, and a decoder module 1490, where each of the modules includes one or more layers of neurons. The encoder module 1410 is configured to receive the microphone signal features {tilde over (D)}_(k,f), the far-end signal features {tilde over (X)}_(k,f), and the error features {tilde over (E)}_(k,f) as input. As noted above, in some embodiments, the recurrent neural network of the echo canceller 1328 merely includes an encoder/decoder GRU network having an encoder module 1410 connected directly to a decoder module 1490 (e.g., without a contextual attention module 1430 in between).

As shown in FIG. 14A, the encoder module 1410 takes the logarithmic magnitudes of the microphone signal features {tilde over (D)}_(k,f)=[D_(t-T+1), . . . , {tilde over (D)}_(t)], the far-end signal features {tilde over (X)}_(k,f)=[{tilde over (X)}_(t-T+1), . . . , {tilde over (X)}_(t)], and the error features {tilde over (E)}_(k,f)=[{tilde over (E)}_(t-T+1), . . . , {tilde over (E)}_(t)] and maps them to hyperspace to generate encoded features h, where the encoded features in the embodiment shown in FIG. 11A have dimensions of T×1000: h=Encoder({tilde over (D)} _(k,f) ,{tilde over (X)} _(k,f) ,{tilde over (E)} _(k,f))

In more detail, in some embodiments of the present disclosure as shown in FIG. 14B, the encoder module 1410 (or encoder layer) includes the concatenation layer 1412 configured to concatenate the input microphone signal features {tilde over (D)}_(k,f), far-end signal features {tilde over (X)}_(k,f), and error features {tilde over (E)}_(k,f) (each having dimensions T×257) to compute concatenated features of dimension T×771. The concatenated features are then supplied to a gated recurrent unit (GRU) layer 1114 substantially similar to the GRU layer shown in FIG. 3 and FIG. 11B, where each unit of the GRU layer 1414 computes the encoded features h in accordance with: h _(t) =f([{tilde over (D)} _(t) ;{tilde over (X)} _(t) ;{tilde over (E)} _(t)],h _(t-1)) where f represents the GRU function and h_(t) represents the output of the GRU for time t.

In embodiments of the present disclosure where the echo canceller 1328 includes a contextual attention neural network, the contextual attention module 1430 takes the encoded features h to compute contextual attentions c that identify important regions of the hyperspace, where the contextual attentions c in the embodiment shown in FIG. 14A have dimensions of T×1000: c=Attention(h)

In some embodiments, the contextual attention module 1430 includes a multi-head self-attention (MHSA) layer and an additional multi-head attention (MHA) layer, substantially similar to that shown in FIG. 11C. Accordingly, the various embodiments of the contextual attention module 1130 described with respect to FIG. 11C may also be used for the contextual attention module 1430 in embodiments of the present disclosure corresponding to a hybrid classical and contextual attention neural network, and the contextual attention module 1430 will not be described in additional detail herein.

Referring to FIGS. 14A and 14B, in the case where the echo canceller 1328 includes a contextual attention neural network, the estimated logarithmic magnitude of the near-end speech {tilde over (Q)}_(k,f) is then computed from the contextual attentions c using the decoder module 1490, where the dimensions of the near-end speech {tilde over (Q)}_(k,f) in the embodiments shown in FIGS. 14A and 14B has dimensions of T×257 (e.g., the same dimensions as the input microphone and far-end signal features). In embodiments of the present disclosure in which the echo canceller 1328 includes an encoder/decoder GRU network, the decoder computes the estimated logarithmic magnitude of the near-end speech {tilde over (Q)}_(k,f) based on the encoded features h.

In more detail, in some embodiments, such as that shown in FIG. 14B, the decoder module 1490 (or decoder layer) includes two GRU layers 1492 and 1494.

According to one embodiment, the first GRU layer 1492 computes its outputs s (s_(k-T+1), s_(k-T+2), . . . , s_(k)) based on the contextual attentions c (c_(k-T+1), c_(k-T+2), . . . , c_(k)) from the contextual attention layer 1430 in accordance with: s _(t) =g ₁(c _(t) ,s _(t-1))

In some embodiments, the first GRU layer 1492 computes its outputs s (s_(k-T+1), s_(k-T+2), . . . , s_(k)) based on the encoded features h (h_(k-T+1), h_(k-T+2), . . . , h_(k)) from the encoder module 1410 in accordance with: s _(t) =g ₁(h _(t) ,s _(t-1)) The second GRU layer 1494 computes the spectral features of the estimated near end speech {tilde over (Q)}_(k,f) ({tilde over (Q)}_(k-T+1); {tilde over (Q)}_(k-T+2), . . . , {tilde over (Q)}_(k)) in accordance with: {tilde over (Q)} _(t) =g ₂(s _(t) ,{tilde over (Q)} _(t-1))

In some embodiments, an exponential linear unit (elu) activation function is used with the first GRU layer 1492 and a linear activation function is used with the second GRU layer 1494.

As shown in FIG. 13, the estimated features of the near-end speech near-end speech {tilde over (Q)}_(k,f) are converted back into estimated near-end speech signal q(t), which may be transmitted to the far-end communication device.

In some embodiments, the encoder/decoder GRU network or the contextual attention neural network 1400 shown in FIGS. 14A, 14B, and 14C is trained using substantially the same training procedure as that described above with respect to FIGS. 5A, 5B, and 5C, and the resulting trained contextual attention neural network 1400 may be deployed into a communication device (e.g., an end-user device such as a smartphone or a tablet computer) in a manner similar to that described above.

FIG. 15 is a flowchart depicting a method for computing an estimated near-end signal from a received far-end signal and a microphone signal using a hybrid of a classical acoustic echo cancellation technique and a deep recurrent neural network in inference mode, according to one embodiment of the present disclosure. As noted above, in some embodiments, the deep recurrent neural network is an encoder/decoder GRU network and in some embodiments the deep recurrent neural network is a contextual attention neural network. As shown in FIG. 15, at 1502, the acoustic echo cancellation system 1300 receives the far-end signal x(t) and, at 1504, the far-end signal feature extraction module 1320 extracts far-end features {tilde over (X)} from the far-end signal x(t).

Similarly, at 1506, the acoustic echo cancellation system 1300 receives the microphone signal d(t) and, at 1508, the microphone near-end signal feature extraction module 1310 extracts microphone signal features {tilde over (D)} from the microphone signal d(t).

At 1510, the acoustic path parameters G are computed using an appropriate classical automatic echo cancellation technique, such as FDNLMS, as discussed above. At 1512, the far-end features {tilde over (X)}, the microphone signal features {tilde over (D)}, and the acoustic path parameters G are used to compute error features {tilde over (E)}.

At 1514, the trained recurrent neural network (e.g., a component of the echo canceller 1328) computes estimated near-end features {tilde over (Q)} from the far-end features {tilde over (X)}, the microphone features {tilde over (D)}, and the error features {tilde over (E)}. As shown, for example, in FIG. 14B, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, features from the context of prior T frames (e.g., the six prior frames k−6, k−5, . . . , k−1) in addition to the current frame k are also supplied to the recurrent neural network. In particular, as shown in FIG. 14B, when computing the estimated near-end features {tilde over (Q)}_(k) for the current frame k, far-end features {tilde over (X)} and microphone features {tilde over (D)} from the current frame k and the T−1 prior frames k−T+1, k−T+2, . . . , k−1 (e.g., far-end features {tilde over (X)}_(k-T+1), {tilde over (X)}_(k-T+2), . . . , {tilde over (X)}_(k), microphone features {tilde over (D)}_(k-T+1), {tilde over (D)}_(k-T+2), . . . , {tilde over (D)}_(k), and error features {tilde over (E)}_(k-T+1), {tilde over (E)}_(k-T+2), . . . , {tilde over (E)}_(k)).

At 1516, feature inversion module 1370 of the acoustic echo cancellation system 1300 computes an estimated near-end signal q(t) for the current frame from the estimated near-end features {tilde over (Q)} of the current frame. As noted above, the features (e.g., the far-end signal features {tilde over (X)}, the microphone features {tilde over (D)}, and the estimated near-end features {tilde over (Q)} may be in a feature space or hyperspace such as STFT space (e.g., spectral features or spectral domain). Accordingly, in some embodiments, the feature inversion module 1370 transforms the estimated spectral features {tilde over (Q)} from the feature space to a time domain signal q(t) suitable for playback on a speaker at the far-end communication device where the feature inversion module 1370 may include an exponential module 1372 and an inverse short-time Fourier transform (iSTFT) module 1374. As shown in FIG. 13, the phase ∠D_(k,f) of the microphone signal d(t) may also be used by the inverse short-time Fourier transform (iSTFT) 1074 when computing the estimated near-end signal q(t).

Experimental Results of Echo Cancellation Using Hybrid Classical and Contextual Attention Neural Network Techniques

In a manner similar to that described above, in the following discussion, a frequency domain normalized least mean square (FDNLMS) (see, e.g., C. Faller and J. Chen, “Suppressing acoustic echo in a spectral envelope space,” IEEE Transactions on Acoustic, Speech and Signal Processing, vol. 13, no. 5, pp. 1048-1062, 2005) is used as a comparative example for measuring the impact of the contextual attention neural network. Another comparative example of an Encoder/Decoder GRU (e.g., the encoder module 1410 and the decoder module 1490 without the contextual attention module 1430) was also used for measuring the impact of the contextual attention module 1430.

In particular, training data was generated from the TIMIT dataset as described above, and then modified in accordance with synthetic and real room impulse responses (RIRs). In the experiments described herein, for the NLMS updating rule, μ was set to 0.5 when there was no double-talk and set to a very small value otherwise. The forgetting factor α was set to 0.6.

The hybrid neural AEC according to some embodiments of the present disclosure was evaluated using a linear model of the acoustic path. In this set of experiments, the real measured RIRs in “office”, “meeting room”, “lecture room”, “stairway1”, “stairway2”, “bathroom”, and lecture room” were used for the training of the hybrid neural AEC, and the “corridor” conditions were used for testing in a phone in hand-held position (HHP). Here, the mismatch between the training RIRs and testing RIR is small as the same recording device was used to capture both RIRs. The similarity between the training and testing RIRs based on average normalized cross correlation (NCC) is 0.97.

Table 9 shows average ERLE values and PESQ scores for the classical (FDNLMS) benchmark, an Encoder/Decoder GRU, and a hybrid encoder/decoder according to some embodiments of the present disclosure:

TABLE 9 ERLE values and PESQ scores in a linear model of acoustic path Testing SER (dB) Method 0 3.5 7 ERLE FDNLMS (Classical) 12.16 11.46 10.52 Encoder/Decoder GRU 55.54 59.36 61.72 Hybrid Encoder/Decoder 56.51 60.49 61.39 PESQ Unprocessed 1.86 2.10 2.33 FDNLMS (Classical) 2.43 2.63 2.81 Encoder/Decoder GRU 2.88 3.10 3.28 Hybrid Encoder/Decoder 2.97 3.16 3.33

As shown in Table 9, the Hybrid Encoder/Decoder outperforms both benchmark AEC models, but outperforms the Encoder/Decoder GRU by a smaller margin, likely because the conditions of the training set and testing set are similar in this experiment, and thus the benefit of the use of the FDNLMS predictor as a component of the hybrid encoder/decoder is reduced.

Accordingly, the performance of embodiments of the present disclosure is further compared to the Encoder/Decoder GRU benchmark under conditions in which the training and testing conditions are different. For this, seven synthetic RIRs were generated for training and again tested on data that was created by a real measured “corridor” RIR. The “corridor” environment was set with reverberation time (T60) selected from {0.2, 0.4, 0.6, 0.8, 0.9, 1.0, 1.25} seconds. Accordingly, the training and testing RIRs in this experiment were significantly different, with the average NCC being about 0.58. The comparison results are given in Table 10, below.

TABLE 10 ERLE values and PESQ scores in a linear model of acoustic path when trained on synthetic RIRs Testing SER (dB) Method 0 3.5 7 ERLE Encoder/Decoder GRU 17.13 20.96 29.68 Hybrid Encoder/Decoder 42.66 47.96 52.47 PESQ Encoder/Decoder GRU 2.43 2.68 2.89 Hybrid Encoder/Decoder 2.76 2.92 3.06

As shown in Table 10, above, in these experiments, the hybrid method according to some embodiments of the present disclosure outperforms the “encoder/decoder GRU” method by a large margin. As such, the hybrid encoder/decoder method according to some embodiments of the present disclosure achieves a superior echo reduction without significant near-end distortion compared against the clean near-end speech.

Additional experiments were performed to study the impact of non-linear models of acoustic paths on embodiments of the present disclosure. In this set of experiments, we used the non-linear (distorted) echo signal y_(nl)(t) in generating the microphone signals, therefore the model contains both power amplifier clipping and loudspeaker distortions. The synthetic RIRs were used for training and the “Corridor” RIR was used for testing. As before, the results were compared against classical FDNLMS alone and the non-hybrid Encoder/Decoder GRU. The results presented in Table 11, below, show that embodiments of the present disclosure outperform the two benchmarks:

TABLE 11 ERLE and PESQ scores in non-linear model of acoustic path and real measured “Corridor” RIR Testing SER (dB) Method 0 3.5 7 ERLE FDNLMS (Classical) 8.16 7.92 7.54 Encoder/Decoder GRU 5.73 5.66 5.91 Hybrid Encoder/Decoder 19.08 19.97 21.64 PESQ Unprocessed 1.79 2.03 2.27 FDNLMS (Classical) 2.16 2.39 2.61 Encoder/Decoder GRU 2.12 2.38 2.64 Hybrid Encoder/Decoder 2.74 2.93 3.09

As such, some embodiments of the present disclosure relate to an architecture for AEC using a hybrid of classical AEC techniques (e.g., FDNLMS) to guide an encoder/decoder GRU network (e.g., a contextual attention neural network) by providing adaptive information. The benefit of a hybrid technique in accordance with embodiments of the present disclosure is especially observed where there is a mismatch between the training conditions used to train the encoder/decoder GRU network and the test conditions (e.g., deployment conditions in the field), and embodiments of the present disclosure provide a performance benefit over comparative benchmark techniques.

FIG. 16 is a block diagram of a near-end device, according to one embodiment of the present disclosure. As noted above, in some embodiments the near-end device 10 is a communications device such as a smartphone or a tablet. The near-end device 10 may include a modem 1610 configured to send and receive data, such as voice communications data. The modem may include a signal processor configured to convert the received encoded data into audio data (e.g., pulse coded modulation data) corresponding to the far-end signal x(t). The decoded data is supplied to a processor 1620, which may implement the acoustic echo cancellation system 200. As noted above, in various embodiments of the present disclosure, the processor implementing the acoustic echo cancellation system 200 may be: a general purpose central processing unit; a graphical processing unit (GPU); a field programmable gate array (FPGA); an neural processing unit (NPU) or neural network processor (NNP) (e.g., a processor having an architecture tailored to perform inference using a neural network); or a neuromorphic processor, although embodiments of the present disclosure are not limited thereto. The received far-end signal x(t) may be played through a loudspeaker 12, which emits sound waves into the environment, some of which may be reflected and received by the near-end microphone 14, which also receives a near-end signal s(t). The generated microphone signal d(t) is supplied to the processor 1620 and the acoustic echo cancellation system 200 to compute the estimated near-end signal q(t), which is transmitted to the far-end device via the modem 1610. Generally, a near-end device 10 may include other components not shown in FIG. 10, such as user interface components (e.g., buttons, a touchscreen display, and the like), a battery or other power source, additional general purpose or special purpose processors, transitory and non-transitory memory, and the like.

Accordingly, aspects of embodiments of the present disclosure relate to recurrent neural networks, including deep multitask recurrent neural networks, for acoustic echo cancellation (AEC). As shown in experimental results, embodiments of the present disclosure perform well in both single-talk and double-talk periods. Some aspects of embodiments of the present disclosure relate to end-to-end multitask learning of both the echo and the near-end signal simultaneously, which improves the overall performance of the trained AEC system. In addition, some aspects of embodiments relate to the use of low-latency causal context windows to improve the context-awareness when estimating the near-end signal with the acoustic echoes removed. When compared based on reference datasets, embodiments of the present disclosure can reduce the echo more significantly than comparative techniques and is robust to additive background noise. Further, a hybrid method according to some embodiments of the present disclosure is more robust to the changes in room impulse response (RIR) and can perform well if fine-tuned by augmenting the data simulated with the impulse response of the target device (e.g., the near-end device 10) under use.

As such, aspects of embodiments of the present disclosure relate to echo cancellation or echo suppression using a trained deep recurrent neural network. While the present disclosure has been described in connection with certain exemplary embodiments, it is to be understood that the disclosure is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the scope of the appended claims, and equivalents thereof. 

What is claimed is:
 1. A system for performing echo cancellation, the system comprising: a processor; and memory storing instructions that, when executed by the processor, cause the processor to: receive a far-end signal from a far-end device at a near-end device; record a microphone signal at the near-end device, the microphone signal comprising: a near-end signal; and an echo signal corresponding to the far-end signal; extract a plurality of far-end features from the far-end signal; extract a plurality of microphone features from the microphone signal; compute a plurality of estimated near-end features from the microphone signal by supplying the microphone features and the far-end features to an acoustic echo cancellation module comprising a contextual attention neural network comprising: an encoder comprising a plurality of gated recurrent units configured to encode the microphone features and the far-end features as encoded features; a contextual attention module configured to compute contextual attention features from the encoded features; and a decoder comprising a plurality of gated recurrent units configured to decode the estimated near-end features from the contextual attention features; compute an estimated near-end signal from the estimated near-end features; and transmit the estimated near-end signal to the far-end device.
 2. The system of claim 1, wherein the encoder comprises one or more of layers of gated recurrent units, and the encoded features correspond to outputs of a last layer of the one or more layers of the gated recurrent units of the encoder, wherein the encoded features are supplied to a multi-head self-attention layer of the contextual attention module; and wherein the decoder comprises one or more layers of gated recurrent units, and the estimated near-end features correspond to an output of a last layer of the one or more layers of gated recurrent units of the decoder.
 3. The system of claim 2, wherein the contextual attention module further comprises a multi-head attention layer configured to receive an output of the multi-head self-attention layer.
 4. The system of claim 3, wherein the contextual attention module further comprises a first layer normalization layer between the multi-head self-attention layer and the multi-head attention layer.
 5. The system of claim 4, wherein the contextual attention module further comprises a second layer normalization layer configured to compute the contextual attention features in accordance with outputs from other layers of the contextual attention module.
 6. The system of claim 1, wherein the memory further stores instructions that, when executed by the processor, cause the processor to: compute a plurality of acoustic path parameters in accordance with prior acoustic path parameters, prior far-end features, and prior microphone features; and compute a plurality of error features based on the far-end features, the microphone features, and the acoustic path parameters, and wherein the acoustic echo cancellation module is further configured to compute the plurality of estimated near-end features from the microphone signal based on the error features.
 7. The system of claim 1, wherein the far-end features supplied to the contextual attention neural network comprise a current frame of far-end features and a causal window of a plurality of previous frames of far-end features, and wherein the microphone features supplied to the contextual attention neural network comprise a current frame of microphone features and the causal window of a plurality of previous frames of microphone features.
 8. The system of claim 7, wherein the estimated near-end features comprise a current frame of estimated near-end features and the causal window of a plurality of previous frames of estimated near-end features, and wherein the contextual attention neural network is trained by iteratively training a plurality of parameters configuring contextual attention neural network to minimize differences between the estimated near-end features and a plurality of ground truth near-end features of the plurality of training data.
 9. The system of claim 8, wherein the training data is generated by: loading a corpus of training data comprising recorded utterances from a plurality of different human speakers; selecting a plurality of pairs of human speakers; for each pair of human speakers: concatenating multiple utterances of a first speaker of the pair of human speakers to generate a training far-end signal; transforming the concatenated utterances to simulate an acoustic path to generate a training echo signal; padding an utterance of a second speaker of the pair of human speakers to generate a training near-end signal having a length equal to the training far-end signal; and mixing the training echo signal with the training near-end signal to generate a training microphone signal.
 10. The system of claim 9, wherein the transforming the concatenated utterances comprises convolving the far-end signal with a room impulse response of a simulated room.
 11. The system of claim 10, wherein the transforming the concatenated utterances further comprises: applying hard clipping to the far-end signal to generate a clipped far-end signal; and applying sigmoidal distortion to the clipped far-end signal.
 12. A system for performing echo cancellation, the system comprising: a processor; and memory storing instructions that, when executed by the processor, cause the processor to: receive a far-end signal from a far-end device; record a microphone signal comprising: a near-end signal; and an echo signal corresponding to the far-end signal; extract a plurality of far-end features from the far-end signal; extract a plurality of microphone features from the microphone signal; compute a plurality of acoustic path parameters in accordance with prior acoustic path parameters, prior far-end features, and prior microphone features; compute a plurality of error features based on the far-end features, the microphone features, and the acoustic path parameters; compute a plurality of estimated near-end features from the microphone signal by supplying the microphone features, the far-end features, and the error features to an acoustic echo cancellation module comprising a recurrent neural network comprising: an encoder comprising a plurality of gated recurrent units configured to encode the microphone features and the far-end features as encoded features; and a decoder comprising a plurality of gated recurrent units configured to decode the estimated near-end features based on the encoded features; compute an estimated near-end signal from the estimated near-end features; and transmit the estimated near-end signal to the far-end device.
 13. The system of claim 12, wherein the acoustic echo cancellation module further comprises a contextual attention module configured to compute contextual attention features from the encoded features, and wherein the decoder is configured to compute the estimated near-end features based on the contextual attention features.
 14. The system of claim 12, wherein the far-end features supplied to the acoustic echo cancellation module comprise a current frame of far-end features and a causal window of a plurality of previous frames of far-end features, wherein the microphone features supplied to the acoustic echo cancellation module comprise a current frame of microphone features and the causal window of a plurality of previous frames of microphone features, and wherein the error features supplied to the acoustic echo cancellation module comprise a current frame of error features and the causal window of a plurality of previous frames of error features.
 15. The system of claim 12, wherein the estimated near-end features comprise a current frame of estimated near-end features and a causal window of a plurality of previous frames of estimated near-end features, and wherein the recurrent neural network is trained by iteratively training a plurality of parameters configuring the recurrent neural network to minimize differences between the estimated near-end features and a plurality of ground truth near-end features of the plurality of training data.
 16. The system of claim 15, wherein the training data is generated by: loading a corpus of training data comprising recorded utterances from a plurality of different human speakers; selecting a plurality of pairs of human speakers; for each pair of human speakers: concatenating multiple utterances of a first speaker of the pair of human speakers to generate a training far-end signal; transforming the concatenated utterances to simulate an acoustic path to generate a training echo signal; padding an utterance of a second speaker of the pair of human speakers to generate a training near-end signal having a length equal to the training far-end signal; and mixing the training echo signal with the training near-end signal to generate a training microphone signal.
 17. The system of claim 16, wherein the transforming the concatenated utterances comprises convolving the far-end signal with a room impulse response of a simulated room.
 18. The system of claim 17, wherein the transforming the concatenated utterances comprises: applying hard clipping to the far-end signal to generate a clipped far-end signal; and applying sigmoidal distortion to the clipped far-end signal.
 19. The system of claim 12, wherein the far-end features, the microphone features, and the estimated near-end features comprise log short time Fourier transform features in logarithmic spectral space. 